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Abstract 

A new algorithm for simulating stochastic partial differential equa- 
tions (SPDEs) of evolutionary type, which is in some sense an in- 
finite dimensional analog of Milstein's scheme for finite dimensional 
stochastic ordinary differential equations (SODEs), is introduced and 
analyzed in this article. The Milstein scheme is known to be impres- 
sively efficient for scalar one-dimensional SODEs but only for some 
special multidimensional SODEs due to difficult simulations of iter- 
ated stochastic integrals in the general multidimensional SODE case. 
It is a key observation of this article that, in contrast to what one 
may expect, its infinite dimensional counterpart introduced here is 
very easy to simulate and this, therefore, leads to a break of the com- 
plexity (number of computational operations and random variables 
needed to compute the scheme) in comparison to previously consid- 
ered algorithms for simulating nonlinear SPDEs with multiplicative 
trace class noise. The analysis is supported by numerical results for a 
stochastic heat equation, stochastic reaction diffusion equations and 
a stochastic Burgers equation showing significant computational sav- 
ings. 

*Arnulf Jentzen, Faculty of Mathematics, Bielefeld University, Germany 
(jentzen@math.uni-bielefeld.de) 

Michael Rockner, Faculty of Mathematics, Bielefeld University, Germany 
(roeckner@math.uni-bielefeld.de) and Department of Mathematics and Statistics, Purdue 
University, USA (roeckner@math.purdue.edu) 



1 



1 Introduction 



Stochastic partial differential equations (SPDEs) of evolutionary type are a 
fundamental instrument for modelling all kinds of dynamics with stochastic 
influence in nature or in man-made complex systems. Since explicit solutions 
of such equations are usually not available, it is a very active research topic 
in the last two decades to solve SPDEs approximatively. The key difficulty 
in this research field is the high computational complexity (in comparison 
to solve deterministic partial and stochastic ordinary differential equations 
approximatively) needed to compute an approximation of the solution of a 
SPDE. The high computational complexity is a consequence of the need 
of discretizing the continued time interval, the infinite dimensional state 
space of the SPDE and the infinite dimensional driving noise process. In 
this article a new numerical method for SPDEs with reduced computational 
complexity in comparison to previously considered algorithms for simulat- 
ing nonlinear SPDEs with multiplicative trace class noise is introduced and 
analyzed. This method is in a sense an infinite dimensional analog of Mil- 
stein's scheme for finite dimensional stochastic ordinary differential equations 
(SODEs). The Milstein scheme is known to be impressively efficient for scalar 
one- dimensional SODEs but only for some special multidimensional SODEs 
due to difficult simulations of iterated stochastic integrals in the general mul- 
tidimensional SODE case. In particular, this suggests that there is no hope 
to expect that an infinite dimensional analog of Milstein's scheme can be 
simulated efficiently and even less to expect that such a method yields a 
break of the computational complexity for approximating SPDEs. However, 
the main contribution of this article is to derive an infinite dimensional ana- 
log of Milstein's scheme that can be simulated very easily by exploiting that 
the infinite dimensional SPDE state space is a function space on which the 
noise acts multiplicatively. In contrast to Milstein's scheme for SODEs, an 
additional ingredient of its infinite dimensional counterpart introduced here 
is a mollifying exponential term approximating the dominant linear opera- 
tor of the SPDE. We also mention that in the case of a linear SPDE, the 
splitting-up method as considered by I. Gyongy and N. Krylov in [20J already 
yields a break of the computational complexity and refer to Section 14.31 for 
a detailed comparison of the splitting-up method and our algorithm in this 
article. In the rest of this introductory section we first review Milstein's 
scheme for SODEs, then reconsider a standard method for solving SPDEs 
approximatively and finally, introduce our algorithm for simulating SPDEs. 

Let T E (0,oo) be a real number, let (ft, J 7 , P) be a probability space 
and let w = (w 1 , . . . , w m ) : [0, T] x tt ->■ R m with m e N := {1,2,...} bea 
m-dimensional standard Brownian motion with respect to a normal filtration 
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(J*t)te[o,T]- Moreover, let deN, x eR d and let /i = (/i^ . . . , /x d ) : R d M d 
and a = (crjj) ig/i,...^} : R d — > M dxm be two appropriate smooth and regular 
je{i,...,m} 

functions with globally bounded derivatives (see, e.g., Theorem 10.3.5 in P. E. 
Kloeden and E. Platen [36] for details). The stochastic ordinary differential 
equation 

dX t = fi(X t ) dt + a(X t ) dw t , X = x (1) 

for all t G [0, T] then admits a unique solution. More precisely, there exists an 
up to indistinguishability unique adapted stochastic process X : [0, T] x Q — > 
M. d with continuous sample paths which satisfies 

X t = x + [ fi{X s )ds+ [ a(X s )dw s (2) 
Jo Jo 

= x + n(X s )ds + Y] / a. t (X s )dwi P-a.s. 
Jo i=1 Jo 

for all t G [0, T]. Here cr, : R d ->■ M d is given by a* (a;) = (01^(2;), . . . , a d ^{x)) 
for all x G W 1 and all i G {1, . . . , m}. Milstein's method (see, e.g., (3.3) in 
Section 10.3 in P. E. Kloeden and E. Platen |36j and also G. N. Milstein's 
original article [13] ) applied to the SODE flTJ is then given by J 7 /B(M d )- 
measurable mappings : Q — > M d , n G {0, 1, ... , N}, iVeN, with t/q^ = x 
and 



Vn+l =Vn + T7 ' Mj/n ) + S ^ 



n 

i=l 



^(n+l)T — W»T 



(n + l)T 



rn a / o \ /» - — 1 — - — /* s 

+ EE ^. /„ r « (3) 

i,j=l k=l 17 Iv 17 lv 

P-a.s. for all n G {0, 1, ... , iV-1} and all iV G N. Although Milstein's scheme 
is known to converge significantly faster than many other methods such as the 
Euler-Maruyama scheme, it is only of limited use due to difficult simulations 

(n + l)T 

of the iterated stochastic integrals f nT N J^r dw 3 u dwl for i,j G {1, . . . , m) 

with i / j, n G {0, 1, ... AT — 1} and iV G N in (J3J). In the special situation 
of so called commutative noise (see (3.13) in Section 10.3 in |36j), i.e. 

s ( ^~ cx ^ ■ akj ^ = s ( ) ^ ■ ^ 

k=l ^ fe ' k=l \ k ' 

for all x G R d and alH, j G {1, . . . , m}, the Milstein scheme can be simplified 
and complicated iterated stochastic integrals in can be avoided. More 
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precisely, in the case (Jl]), Milstein's scheme (J3J) reduces to 
T m ( \ 

Vn+l =Vn ' + jy ' Mfn ) + E ) ' ( W '"+ 1 ) T ~ J 

i=l 

+ ^ E E (^) ■ ) • (^w - <%) • (^w - 

i,i=l fc=l 7 

i=l k=l V K 7 

P-a.s. for all n G {0, 1, . . . , JV - 1} and all JV G N (see (3.16) in Section 10.3 
in [US])- For instance, in the case d = m = 1, condition ([!]) is obviously 
fulfilled and the Milstein scheme (JS} can then be written as 

Vn+l = Vn + ]y ■ Ml/n ) + ) ' ~ w ^f) 

+ i • ^(v?) • cr(^) • (("feflir - - |) (6) 

P-a.s. for all n G {0, 1, . . . , JV - 1} and all JV G N (see (3.1) in Section 10.3 
in [36]). Of course, (jUJ) can be simulated very efficiently However, (jl]) is 
in the case of a multidimensional SODE seldom fulfilled and even if it is 
fulfilled, Milstein's method (jSJ becomes less useful if d, m G N are large. 
For example, if d = m = 20 holds, then the middle term in fl5]) contains 
20 3 = 8000 summands. So, more than 8000 additional arithmetic operations 
are needed to compute y^+i 

from y% for n G {0, 1, . . . , JV - 1}, JV G N 
via (|SD in the case d = m = 20 in general which makes Milstein's scheme 
less efficient. This suggests that there is no hope to expect that an infinite 
dimensional analog of Milstein's method can be simulated efficiently in the 
case of infinite dimensional state spaces such as L 2 ((0, 1),R) instead of M. d 
and M m respectively. The purpose of this article is to demonstrate that 
this is not true. More precisely, an infinite dimensional analog of Milstein's 
scheme that can be simulated very easily is derived and analyzed here. This 
leads to a break of the complexity (number of computational operations and 
random variables needed to compute the scheme) in comparison to previously 
considered algorithms for simulating nonlinear SPDEs with multiplicative 
trace class noise which is illustrated in the following. 

Let H = L 2 ((0,1),M) be the M-Hilbert space of equivalence classes of 
£>((0, l))/i3(M)-measurable and Lebesgue square integrable functions from 
(0, 1) to R and let /, b : (0, 1) x M. — > M. be two appropriate smooth and 
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regular functions with globally bounded derivatives (see f[3"0"j) and (jHTl) for 
details). As usual we do not distinguish between a B((0, 1))/ i3(lR)-measurable 
and Lebesgue square integrable function from (0, 1) to R and its equivalence 
class in H . Moreover, let k G (0, oo) be a real number, let £ : [0, 1] — > R 
with £(0) = = be a smooth function and let W : [0,T] x Q -> 

if be a standard Q- Wiener process with respect to F t , t G [0,T], with a 
trace class operator Q : H ^ H (see, for instance, Definition 2.1.12 in 
[50J). It is a classical result (see, e.g., Proposition 2.1.5 in (50]) that the 
covariance operator Q : H — > H of the Wiener process W : [0, T] x — )■ H 
has an orthonormal basis Qj £ H, j £ N, of eigenfunctions with summable 
eigenvalues /ij G [0, oo), j G N. In order to have a more concrete example, 
we consider the choice gj(x) = y/2sm(jnx) and fxj = j% for all x G (0, 1) and 
all j G N in the following and refer to Section [2] for our general setting and 
to Section H] for further possible examples. Then we consider the SPDE 



dXt(x) 



K -^X t (x) + f(x,X t (x)) 



dt + b(x,X t (x))dW t (x) (7) 



with X t {0) = X t {l) = and X (x) = £{x) for x G (0,1) and t G [0,T] 
on H. Under the assumptions above the SPDE (j7j) has a unique mild solu- 
tion. Specifically, there exists an up to indistinguishability unique adapted 
stochastic process X : [0, T] x $7 — >• H with continuous sample path which 
satisfies 

X t = e At i+ [ e A{t - s) F(X s )ds+ [ e A{t - s) B(X S ) dW s P-a.s. (8) 
Jo Jo 

for all t G [0, T] where A : -D(A) C H — > H is the Laplacian with Dirichlet 
boundary conditions times the constant k G (0, oo) and where F : H — > 
H and B : H — >■ HS(U ,H) are given by (F(t>))(x) = f(x,v(x)) and 
(-B(f)-u)(x) = 6(x,f(x)) • m(x) for all x G (0,1), v £ H and all w G C/o- 

Here ?7o = Q^{H) with (v,w) Uq = (q~^v,Q~^w^ for all v,w & Uo is the 

image M-Hilbert space of (see Appendix C in [HO] ) . 

Then our goal is to solve the strong approximation problem (see Sec- 
tion 9.3 in [36J) of the SPDE (J7|). More precisely, we want to compute a 
F / B(H)-measmab\e numerical approximation Y : Q — > H such that 



E 



i 

\X T (x) — Y(x)\ 2 dx 



2 



<e (9) 

L^o 

holds for a given precision e > with the least possible computational effort 
(number of computational operations and independent standard normal ran- 
dom variables needed to compute Y : Q — > H). A computational operation is 
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here an arithmetic operation (addition, subtraction, multiplication, division), 
a trigonometric operation (sine, cosine) or an evaluation of / : (0, 1) xR — >■ R, 
b : (0, 1) x R — > R or the exponential function. 

In order to be able to simulate such a numerical approximation on a 
computer both the time interval [0, T] and the infinite dimensional space 
H = L 2 ((0, 1),R) have to be discretized. While for temporal discretizations 
the linear implicit Euler scheme (see |3 El UHl ESI I2S1 EZl EZl ESI EO]) and 
the linear implicit Crank-Nicolson scheme (see [25j |26j EH [58]) are often 
used, spatial discretizations are usually achieved with finite elements (see 
[H El El El HH1 EZl EH EE1 [391 1121 EH1 ED]), finite differences (see [HH E2 H31 
US E21 [531 [52 EH E3 [59]) and spectral Galerkin methods (see [151 ESI EHl ED 
EZl HH El EH SB])- For instance, the linear implicit Euler scheme combined 
with spectral Galerkin methods which we denote by J-"/i3(if)-measurable 
mappings Z* : O ->• H, n G {0, 1, . . . , N 3 }, N G N, is given by = P N (£) 
and 

=Pn(i~ ^ + ^ • /(., Z2) + b(; Z») ■ (W&pr ~ W^J) 

T (10) 

P-a.s. for all n G {0, 1, ... , A^ 3 — 1} and all N EN. Here the bounded linear 
operators Pn : H -+ H, N e N, and the Wiener processes W N : [0, T] x Q — > 
H , A" G N, are given by 

N „i 

(Pn(v))(x) — 2 sin(n7ra;) / sm(niiy) v(y) dy 

n=l ^ 

for all x G (0,1), v G if, JV G N and by W^w) = *V(W t (a;)) for all 
t G [0, T], u E Q, N E N. Moreover, we use the notations v • w : (0, 1) — > R, 
v 2 : (0, 1) -»■ R and y>(-, u) : (0, 1) R given by 

(y ■ w)(x) = v(x) ■ w(x), (w 2 )(x) = (v(x)) 2 , (</?(•, v))(x) = <p(x, v(x)) 

for all x G (0, 1) and all functions v, w : (0, 1) — > R, if : (0, 1) x R — > R here 
and below. In ([TO]) the infinite dimensional R-Hilbert space H is projected 
down to the N- dimensional R-Hilbert space Pn(H) with iV G N and the 
infinite dimensional Wiener process W : [0, T] X O — > H is approximated 
by the finite dimensional processes W N : [0, T] x fi — > H, N G N, for 
the spatial discretization. For the temporal discretization in the scheme 
Z^ , n G {0, 1, . . . , iV 3 }, above the time interval [0,T] is divided into iV 3 
subintervals, i.e. A^ 3 time steps are used, for A" G N. The exact solution 
X : [0,T] x Q -»■ H of the SPDE © has values in D((-A)t) and satisfies 
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E\\(-A)~<X T \\ 2 H < oo for all 7 G (0, f) (see Section 4.3 in [32]). This shows 
(E \\X T - Pn{X t )\\ 2 h Y < (E \\{-Ayx T \\ 2 H y \\(-A)^(I - P N )\ 



\L(H) 



< (E ||(-A) 7 X r ||y 5 (1 + fiT 1 ) A^ -27 < oo 



for all N G N and all 7 G (0, §). So, P N {X T ) converges in the root mean 
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square sense to X? with order |— as A/" goes to infinity. (For a real number 
5 G (0,oo), we write 5— for the convergence order if the convergence order 
is higher than 8 — r for every arbitrarily small r G (0, 8).) Additionally, the 
solution process of the SPDE (J7]) is known to be —Holder continuous in the 
root mean square sense (see, for instance, Theorem 1 in [34]) and therefore, 
the linear implicit Euler scheme converges temporally in the root mean square 
to the exact solution of the SPDE (j7j) with order ~ (see, e.g., Theorem 1.1 
in [60]). Combining the convergence rate |— for the spatial discretization 
and the convergence rate | for the temporal discretization indicates that it 
is asymptotically optimal to use the cubic number A^ 3 of time steps in the 
linear implicit Euler scheme , ue {0,1,..., A^ 3 }, above. 

We now review how efficiently the numerical method ffTUj) solves the strong 
approximation problem d^J) of the SPDE ([7j). Standard results in the litera- 
ture (see, for instance, Theorem 2.1 in [26]) yield the existence of real numbers 
C r > 0, r G (0, |), such that 



E 



\X T (x) — Z^ 3 (x) 1 2 dx 



<C r -N (r ~^ (11) 



holds for all A^ G N and all arbitrarily small r G (0, |). The linear implicit 
Euler approximation Z^ 3 thus converges in the root mean square sense to X? 
with order |— as A" goes to infinity. Moreover, since Pn(H) is A^-dimensional 
and since A^ 3 time steps are used in ffTOl . 0(N 4 log(AT)) computational opera- 
tions and random variables are needed to compute Z^ :i . The logarithmic term 
in 0(N 4 \og(N)) arises due to computing the nonlinearities / and b with fast 
Fourier transform where aliasing errors are neglected here and below. Com- 
bining the computational effort 0(N A log(A^)) and the convergence order | — 
in pip shows that the linear implicit Euler scheme needs about 0(e _ §) com- 
putational operations and independent standard normal random variables to 
achieve the desired precision e > in Q. In fact, we have demonstrated that 
Euler's method (TlOjl needs 0(e _( -3 +T ')) computational operations and random 
variables to solve ([9]) for every arbitrarily small r G (0, 00) but for simplicity 
we write about 0(e~z) computational operations and random variables here 
and below. 
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Having reviewed Euler's method flTUj) . we now derive our infinite dimen- 
sional analog of Milstein's scheme. In the finite dimensional SODE case, the 
Milstein scheme fl3]) is derived by applying ito's formula to the integrand 
process a(X t ), t G [0, T], in t^j. This approach is based on the fact that the 
diffusion coefficient a is a smooth test function and that the solution process 
of ([T]) is a Ito process. This strategy is not directly available in infinite di- 
mensions since (|7|) does in general not admit a strong solution to which Ito's 
formula could be applied. Recently, in [32] in the case of additive noise and 
in [30] in the general case, this problem has been overcome by first applying 
Taylor's formula in Banach spaces to the diffusion coefficient B in the mild 
integral equation (jSJ) and by then inserting a lower order approximation re- 
cursively (see Section 4.3 in [3U]). More formally, using F(X S ) ~ F(X ) and 
B(X S ) « B(X ) + B'(X )(X S - X ) for s G [0,T] in © shows 

X t « e At i + f e A{t - s) F{X ) ds + / e A{t - s) B{X ) dW s 
Jo Jo 

+ [ e A ^B'(X )(X s -X )dW s 
Jo 

« e At (x + t- F(X ) + J B(X ) dW s + J B'(X )(X S - X ) dw)j 

for t G [0, T). The estimate X s « X + j° B(X ) dW u for s G [0, T] then gives 

X t « e At (jf + 1 ■ F(X ) + jf B(X ) dW s + jf 5'(X ) Qf 5(X ) tffl^ cW^ 

(12) 

for t G [0, T]. Using Ito's formula this temporal approximation has already 
been obtained in (1.12) in [IS] under additional smoothness assumptions 
of the driving noise process of the SPDE (jHJ) (see Assumption C in |45j) 
which guarantee the existence of a strong solution and thus allow the appli- 
cation of Ito's formula. Combining the temporal approximation Q12j) and the 
spatial discretization in ffTU]) indicates the numerical scheme with J- '/ 13(H)- 
measurable mappings : Q — > H, n G {0, 1, ... , N 2 }, N e N, given by 
Y N = P N (0 and 

*£i = ^ + ^ • F(^) + 5(1^) - W|) 

^iV;; N )( / IMY^dwAdW") (13) 
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P-a.s. for all n G {0, 1, . . . , iV 2 — 1} and all N G N. Now, we are at a stage 
similar to the finite dimensional case a higher order method seems to 
be derived which nevertheless seems to be of limited use due to the iterated 
high dimensional stochastic integral in ( fl3|) . However, a key observation of 
this article is the following formula 

(n + l)T / \ 

P* - B'(Y n N ) f B(Y n N ) dW u N dW» 
= \ {ly b ) ^ • ^ • ((*% ~ W i) 2 - ^ £ 

(14) 

P-a.s. for all n G {0, 1, . . . , iV 2 — 1} and all JV G N (see Subsection O for 
the proof of the iterated integral identity (1T41) and see below for a heuristic 
explanation of this fact). So, the iterated high dimensional stochastic integral 
in f JT3|l reduces to a simple product of functions. The function J^b : (0, 1) x 
R — ?• R is here the partial derivative (J^b)(x,y) for x G (0, 1) and y G R. 
Using (fT4"]) the numerical scheme ( IT5|) thus reduces to 

= Pn e A ^ (r n N + ^ • /(, O + b(; Y n N ) ■ (w&gr - W^j 

+ \ ^ ■ ^ •(("%- w l) 2 - ) 

(15) 

P-a.s. for all n G {0, 1, . . . , iV 2 — 1} and all iV G N. It can be seen that only 
increments of the finite dimensional Wiener processes W N : [0, T] x Q — > H, 
N G N, are used in (fl5|) . Note that, as in the case of ffTOj) . the infinite 
dimensional R-Hilbert space H is projected down to the N- dimensional R- 
Hilbert space Pn(H) with iVeN and the infinite dimensional Wiener process 
W : [0, T] x — y H is approximated by the finite dimensional Wiener pro- 
cesses W N : [0, T] x D, ->■ if, iV G N, for the spatial discretization in ffTol) . 
For the temporal discretization in the scheme , n G {0, 1, ... , A^ 2 }, above 
the time interval [0,T] is divided into iV 2 subintervals, i.e. iV 2 instead of A^ 3 
time steps are used in (Tl5|) . for A^ G N. In the following we explain why it is 
crucial to use A^ 2 time steps in (fl~5"j) instead of A^ 3 time steps in the case of 
the linear implicit Euler scheme fllOp . 

More formally, we now illustrate how efficiently the method ffT5l) solves 
the strong approximation problem (Q of the SPDE ([7]). Theorem [1] (see 
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Section [3] below) gives the existence of real numbers C r > 0, r G (0, |), such 
that 

i 



E 



i 

— Y^>(x)| cfx 



L^o 



3 



< C r ■ (16) 



holds for all N G N and all arbitrarily small r G (0, |). The approximation 
Y$2 thus converges in the root mean square sense to Xt with order |— as iV 
goes to infinity. The expression 

s (IO ( '' FjV) ' K,YnN) ' " ^ " w^ i{9i)2 ) (17) 



8=1 



for n G {0,1,..., X 2 — 1} and N G N in (fl5l) contains additional informa- 
tion of the solution process of ((7|) and this allows us to use less time steps, 
N 2 in (I15p instead of N 3 in (TT0]) . to achieve the same convergence rate as 
the linear implicit Euler scheme (ITU]) (compare (Hip and (H6p ). Nonethe- 
less, ( II 7p and hence the numerical method (1151) can be simulated very eas- 
ily. The function ™ ^2i=i fa (di) 2 i n dlZJ) can be computed once in ad- 
vance for which 0(N 2 ) computational operations are needed. Having com- 
puted -^2 J2iLi fa (.9i) 2 > 0(Nlog(N)) further computational operations and 
random variables are needed to compute (|17|) from for one fixed n G 
{0, 1, . . . , N 2 — 1} by using fast Fourier transform. Since 0(Nlog(N)) com- 
putational operations and random variables are needed for one time step and 
since N 2 time steps are used in (Tl5i) . 0(N 3 log(iV)) computational operations 
and random variables are needed to compute Y^ 2 • Combining the computa- 
tional effort 0(N 3 log(iV)) and the convergence order |— in ( Fl6|) shows that 
the numerical method (j!5p needs about 0(e~ 2 ) computational operations 
and independent standard normal random variables to achieve the desired 
precision e > in (jHJ). To sum up, the algorithm (fT5l) reduces the complexity 
of the problem © of the SPDE from about 0(e _ f ) to about 0(e" 2 ). 

However, the complexity rates 0(e~z) and 0(s~ 2 ) are both asymptotic 
results as e > tends to zero. Therefore, from a practical point of view, 
one may ask whether the algorithm ( 1151) solves the strong approximation 
problem more efficiently than the linear implicit Euler scheme ( II Op for a 
given concrete e > and a given example of the form (J7|) . In order to analyze 
this question we compare both methods in the case of a simple stochastic 
reaction diffusion equation. More formally, let K — let £ : [0, 1] — y R be 
given by £(x) = for all x G [0, 1] and suppose that /, b : (0, 1) x R — > R are 
given by f(x,y) = 1 — y and b(x,y) = for all x G (0, 1), y G R. The 
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SPDE ([7D thus reduces to 




t (x) + 1 -X t (x) dt + 



l-X t {x) 

i + x t { x y 



dW t (x) 



(18) 



with X t (0) = X t (l) = and X = for x G (0, 1) and t G [0, 1] (see also 
Section I4TT1 for more details concerning this example). Additionally, assume 
that ([9]) for the SPDE (j!8p should be solved with the precision of say three 
decimals, i.e. with the precision e = in fl9]). In Figure Q] the approxi- 
mation error in the sense of fl9]) of the linear implicit Euler approximation 
Z^ 3 (see fTTO]) ) and of the approximation Y^ 2 (see ([15]) ) is plotted against the 
precise number of independent standard normal random variables needed to 
compute the corresponding approximation for N G {2,4,8,16,32,64,128}: 
It turns out that in the case of the linear implicit Euler scheme f lTU]) 

and that * n ^ ne case °f ^he algorithm (fl5|) achieve the desired precision 
s = ^ j n f or the SPDE (USD- The Matlab codes for simulating Z™ 3 
via ( ITU]) and K^gl via ( !T5|) for the SPDE (1X8]) are presented below in Figure [2] 
and Figure [3] respectively. The differences of the codes and the additional 
code needed for the new algorithm ( 115)) are printed bold in Figure El The 
Matlab code in Figure |2] requires on an Intel Pentium D a CPU time 
of about 15 minutes and 25.03 seconds (925.03 seconds) while the code 
in Figure [3] requires a CPU time of about 8.93 seconds to be evaluated on 
the same computer. So, the algorithm (fT5"|) is for the SPDE ( TTgj) more than 
hundred times faster than the linear implicit Euler scheme ( ITU]) in order 
to achieve a precision of three decimals in (J9j). To sum up, the four code 
changes (line 1, lines 3-6, line 10 and line 11) in Figure [3] enables 
us to simulate ( 118]) with the precision of three decimals in about 9 
seconds instead of about 15 minutes. Further numerical examples for 
the algorithm ( TT5"j) can be found in Section 0] and Section 

Having illustrated the efficiency of the method f lT5|) . we now take a short 
look into the literature of numerical analysis for SPDEs. First, it should be 
mentioned that any combination of finite elements, finite differences or spec- 
tral Galerkin methods for the spatial discretization and the linear implicit Eu- 
ler scheme or also the linear implicit Crank-Nicolson scheme for the temporal 
discretization do not reduce the complexity 0(e~») of the problem ([9D of the 
SPDE (J7J). However, in the case of a linear SPDE, the splitting- up method as 
considered by I. Gyongy and N. Krylov in [20] (see also H [T2j UM EEJ E2J [23] 
and the references therein) converges with a higher temporal order and there- 
fore breaks the computational complexity in comparison to the linear implicit 
Euler scheme. The key idea of the splitting-up method is to split the consid- 
ered SPDE into appropriate subequations that are easier to solve than the 
original SPDE, e.g., that can be solved explicitly. The splitting-up method 



11 




Precise number of used random variables 



Figure 1: SPDE f[T8"j) : Approximation error in the sense of OH]) of the lin- 
ear implicit Euler approximation (see (ITU]) ) and of the approximation 
Y^ 2 (see (fl5|) ) against the precise number of independent standard normal 
random variables needed to compute the corresponding approximation for 
N E {2,4,8,16,32,64,128}. 

thus essentially depends on the simplicity of the splitted subequations and 
can therefore in general not be used efficiently for nonlinear SPDEs which 
are investigated in this article. Nonetheless, in the case of a linear SPDE as 
in I. Gyongy and N. Krylov's article |20j, the splitting- up method and the 
scheme in this article converge with the same complexity rate (see Section [4731 
for a more detailed comparison of the splitting-up method and the algorithm 
in this article). Additionally, in the case / = in (J7J), T. Miiller-Gronbach 
and K. Ritter invented a new scheme which reduces the number of random 
variables needed for solving a similar problem as ([9]) from about 0(e~^) to 
about 0(e~ 2 ) (see [46J and also [IZ]). Nonetheless, the number of compu- 
tational operations needed and thus the overall computational complexity 
could not be reduced by their algorithm. Moreover, Milstein type schemes 
for SPDEs have been considered in [5] [151 S3 US]. In [IS], W. Grecksch 
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N= 128; M=N~3; A = -pi " 2 * ( 1 :N) . " 2 / 1 00 ; Y= zeros(l,N); 
mu = (l:N).*-2; f = @(x) 1-x; b = @(x) (1-x) ./(1+x . " 2 ) ; 
for m= 1 M 

y = dst (Y) * sqrt (2) ; 

dW = dst( randn(l,N) .* sqrt (mu*2/M) ); 

y = y + f(y)/M + b( y ).*dw ; 

Y = idst ( y ) / sqrt (2) ./ ( 1 - A/M ) ; 
end 

plot( (0:N+1)/(N+1), [0 ,dst(Y)*sqrt (2) ,0] ); 

Figure 2: Matlab code for simulating the linear implicit Euler approxima- 
tion Z% 3 with N = 128 (see (HQJ) for the SPDE JTgJ. 

and P. E. Kloeden proposed a Milstein like scheme for a SPDE driven by a 
scalar one-dimensional Brownian motion (see also [37]). In view of ([6]), their 
Milstein scheme can be simulated efficiently since the driving noise process 
is one-dimensional. Furthermore, in the case of a linear SPDE, P.-L. Chow 
et al. constructed in the interesting article [5] a scheme similar to ffTB"]) but 
with an additional term. (The additional term may be useful for decreasing 
the error constant but turns out not to be needed in order to achieve the 
higher approximation order due to Theorem Q] here.) In order to simulate 
the iterated stochastic integral in their scheme, they then suggest to omit 
the summands in the double sum in (JSJ) for which i ^ j holds (see (2.8) in 
[5]). Their idea thus yields a scheme that can be simulated very efficiently 
but does in general not converge with a higher order anymore except for a 
linear SPDE driven by a scalar one-dimensional Brownian motion. Finally, 
based on Ito's formula, Y. S. Mishura and G. M. Shevchenko proposed in [15] 
the temporal approximation (fT2l under additonal smoothness assumptions 
of the driving noise process of the SPDE (|SJ) which garantuee the existence 
of a strong solution and thus allow the application of Ito's formula (see As- 
sumption C in |45j). The simulation of the iterated stochastic integrals in 
their numerical approximation remained an open question (see Remark 1.1 in 
[35]). To sum up, to the best of our knowledge the computational complexity 
barrier 0(e~») for the problem ([9]) of the SPDE (J7J has not been broken by 
any algorithm proposed in the literature yet. 

It is in some sense amazing that the Milstein type scheme f|T3|) reduces 
to the simple algorithm ffTB]) in the infinite dimensional SPDE setting (jTJ 
which is the key observation in order to beat the computational complexity 
barrier 0(e~a). The reason for this simplification is that B(v) G HS(Uq, H) 
for v G H does not act as an arbitrary linear operator on the infinite di- 
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N= 128; M = N~2; A = -pi * 2 * ( 1 :N) . " 2 / 1 00 ; Y= zeros(l,N); 
mu = (l:N).~-2; f = @(x) 1-x; b = @(x) (l-x)./(l + x/2); 
bb = @(x) (l-x).*(x.~2-2*x-l)/2./(l+x.~2).*3; g = zeros ( 1 ,N) ; 
for n—lN 

g = g+2*sin(n*(l:N)/(N+l)*pi).~2*mu(n)/M; 
end 

for m= 1 M 

y = dst (Y) * sqrt (2) ; 

dW = dst( randn(l,N) .* sqrt (mu*2/M) ); 

y = y + f(y)/M + b( y ).*dw + bb( y ).*(dw.~2 - g ) ; 

Y = exp( A/M ) .* id st ( y ) / sqrt (2) ; 
end 

plot( (0:N+1)/(N+1), [0,dst(Y)*sqrt (2) ,0] ); 

Figure 3: Matlab code for simulating the approximation Y^ 2 with N = 128 
(see (USD) for the SPDE (Jig). 

mensional vector space Uq but as a multiplication operator on the function 
space U Q C L 2 ((0, 1),M). First, this ensures that B : H -> HS(U ,H) 
naturally falls into the infinite dimensional analog of the commutative noise 
case © although b : (0, 1) x E ->■ M and hence B : H ->■ HS(U , H) are 
possibly nonlinear mappings (see AM]) for details). Second, in the infinite 
dimensional commutative noise setting, the acting of B(v) for v G H as a 
multiplication operator assures that the infinite dimensional analog of §J5§ 
even simplifies in infinite dimensions to ffT5l) . These two facts guarantee that 
the method f|T3l) reduces to the simple algorithm (ITS]) which is in a sense the 
main contribution of this article. To sum up, due to the reduced complexity, 
the algorithm (ITS]) provides an in comparison to previously considered nu- 
merical methods impressive instrument for simulating nonlinear stochastic 
partial differential equations with multiplicative trace class noise. 

The rest of this article is organized as follows. In Section [2]the setting and 
the assumptions used are formulated. The numerical method and its con- 
vergence result (Theorem |TJ) are presented in Section [3J In Section H] several 
examples of Theorem UJ including a stochastic heat equation and stochastic 
reaction diffusion equations are considered. Although our setting in Section [2] 
uses the standard global Lipschitz assumptions on the nonlinear coefficients 
of the SPDE, we demonstrate the efficiency of our method numerically for 
a stochastic Burgers equation with a non-globally Lipschitz nonlinearity in 
Section [SJ The proof of Theorem UJ is postponed to Section [HI 

Finally, we would like to add some concluding remarks. There are a num- 
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ber of directions for further research arising from this work. One is to analyze 
whether the exponential term in (|15[) can be replaced by a simpler mollifier 
such as (J — jpA)^ 1 for N G N. This would make the scheme even simpler 
to simulate. A second direction is to combine the temporal approximation in 
(fl5|) with other spatial discretizations such as finite elements. This makes it 
possible to handle more complicated multidimensional domains on which the 
eigenf unctions of the Laplacian are not known explicitly. A third direction 
is to reduce the Lipschitz assumptions in Section [2] in order to handle the 
stochastic Burgers equation in Section[5]and further SPDEs with non-globally 
Lipschitz nonlinearities such as the stochastic porous medium equation and 
hyperbolic SPDEs. Finally, a combination of Giles' multilevel Monte Carlo 
approach in [JJ] with the method in this article should yield a break of the 
computational complexity in comparison to previously considered algorithms 
for solving the weak approximation problem (see, e.g., Section 9.4 in [36]) of 
the SPDE ©. 

2 Setting and assumptions 

Throughout this article suppose that the following setting and the follow- 
ing assumptions are fulfilled. Fix T G (0, oo), let P) be a proba- 
bility space with a normal filtration (J : t)t<=\o,T\ and let (H, (•, ■) H , \\-\\ H ) and 
(U, (•, ■) u , W'Wu) be two separable R-Hilbert spaces. Moreover, let Q : U — > U 
be a trace class operator and let W : [0, T] x Q — > U be a standard Q- Wiener 
process with respect to {J r t)te[o,T]- 

Assumption 1 (Linear operator A). Let I be a finite or countable set and 
let (\) ieX C (0, oo) be a family of real numbers with infj 6 jAj G (0, oo). 
Moreover, let (ei)i e x be an orthonormal basis of H and let A : D(A) C H — > 
H be a linear operator with 

Av = (ei,v) H ei (19) 

for every v G D(A) and with D(A) = {w G H\ J2iei l^*| 2 I ( e *' w )h\ 2 < °°} ■ 

By V r := D ((— A) r ) equipped with the norm \\v\\ v ^ := || (— A) r v\\ H for 
all v G V r and all r G [0, oo) we denote the P-Hilbert spaces of domains of 
fractional powers of the linear operator —A : D(A) C H — >■ H . 

Assumption 2 (Drift term F). Let (3 G [0,1) be a real number and let 
F : Vp — >■ H be a twice continuously Frechet differentiable mapping with 

sup„ e va \\ f '( v )\\l(h) < 00 and su Pt-e^ \\ f "( v )\\lW(Vp,h) < °°- 
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In order to formulate the assumption on the diffusion coefficient of our 
SPDE, we denote by (U , (-, -) Uo , ||-||{/ ) the separable M-Hilbert space U : = 

Q^{U) with (v,w) Uo = (^Q~^v,Q~^w^ for all v,w G U$ (see, for example, 

Section 2.3.2 in |50j). For an arbitrary bounded linear operator S G L(U), we 
denote by S" 1 : im(S') C U — > U the pseudo inverse of S (see, for instance, 
Appendix C in [50J). 

Assumption 3 (Diffusion term B). Let B : Vp — > HS(Uq, H) be a twice con- 
tinuously Frechet differentiable mapping with sup veV \\B'( v )\\l(h hs(u h)) < 
oo and sup veV0 \\B"(v)\\ L(2){v ^ HS{UotH)) < oo. Moreover, let a, c G (0,oo), 

5,$ G (0, |), 7 G [max(5, 0) , 8 + |) be real numbers, let B(Vs) C HS(U ,V S ) 
and suppose that 

\\B(u)\\ HS(Uoys) <c(l+\\u\\ Vs ), (20) 
\\B'(v) B(v) - B'(w) B(w)\\ HS(2HUojH) <c\\v-w\\ H , (21) 

\{-A)-*B(v)QT 

holds for every u & Vs and every v,w G V 1 . Finally, let the bilinear Hilbert- 
Schmidt operator B'{v)B{v) G HS^ 2 \U , H) be symmetric for every v G Vp. 

The operator B'(v)B(v) : Uq x Uq — > H given by 

(B'(v)B(v)) (u,u) = {B'{v) (B(v)u)} (u) 

for all u,u G Uq is a bilinear Hilbert-Schmidt operator in HS^(Uq, H) = 
HS(Uo eg) Uo, H) for every v G Vg. Therefore, condition ( 12T|) means that the 
mapping 

V 1 HS {2) (U , H), v i — y B'(v)B(v), v G V 7 , 
is globally Lipschitz continuous with respect to and ||'|lHS'( 2 )(f7o h) % 

We 

also would like to point out that the assumed symmetry of B'(v)B(v) G 
HS^ 2 \Uq,H) for all v G Vp is the abstract (infinite dimensional) analog of 
flU). More formally, if H = R d , U = R m and Q = I with d, m G N holds, then 
the symmetry of B'{v)B{v) G HS^iR" 1 , R d ) reduces to (gj) (with a replaced 
by B). Although (J3J) is seldom fulfilled for finite dimensional SODEs, the 
symmetry of B'{v)B{v) G HS^^Uq, H) for all v G Vp is naturally met in the 
case of multiplicative noise on infinite dimensional function spaces as we will 
demonstrate in Sections H] and El Finally, we emphasize that we assumed in 
no way that the linear operators A : D(A) C H — )■ H and Q : U — >■ U are 
simultaneously diagonalizable. 



Hsa/n.m 



< c 1 + llf I 



(22) 



16 



Assumption 4 (Initial value £). Let £ : Q — > be a J-q/B (V~) -measurable 



7 u ^ v •> 0/ V y J J 

< oo. 

7 

These assumptions suffice to ensure the existence of a unique solution of 



mapping with E ||£||^ < oo 



the SPDE (im 



Proposition 1 (Existence of the solution). Let Assumptions U§]\ in Sec- 
tion^ be fulfilled. Then there exists an up to modifications unique predictable 
stochastic process X : [0,T] x Q Vy which fulfills sup tg [ 0T ] E ||^t||y 7 < oo ; 

su Pte[o,T] E \\ B ( x t)\\% S (u Q ,v s ) < 00 and 



X t = e At £ + 



f e A{t - s) F(X s )ds+ [ e A{t - s) B(X s )dW s F-a.s. (23) 
Jo Jo 

for all t G [0,T]. Moreover, we have 

sup 1 " t2 . <oo (24) 

hy£t 2 

for every r G [0, 7]. 

Proposition [T] immediately follows from Theorem 1 in 



3 Numerical scheme and main result 

In this section our numerical method is introduced and its convergence result 
is stated. To this end let J be a finite or countable set, let (9j)j e j C U be an 
orthonormal basis of eigenfunctions of Q : U —tU and let C [0, 00) 

be the corresponding family of eigenvalues (such an orthonormal basis of 
eigenfunctions exists since Q : U — > U is a trace class operator, see Proposi- 
tion 2.1.5 in [50J). In particular, we have 

Qu = ^2nj (g j ,u) u g j (25) 

for all u G U. Additionally, let (ZN) NeN and {Jk)k&, De sequences of finite 
subsets of X and J respectively. Then we define the linear projection oper- 
ators P/v : H — > H, N G N, by Pn{v) '■— J2iei N { e ii v )u e i f° r an v e H anc ^ 
all iV G N. Furthermore, we define Wiener processes W K : [0, T] x O — > U$, 
K G N, by 

W t %;):= J2 lSi>WM)v9j (26) 
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for all t £ [0,T], wen and all K £ N. We also use the J r /S([/ )-nieasurable 
mappings AW^' K : -+ E7 , m £ {0, 1, . . . , M — 1}, M,K £ N, given by 
AW%> k (oj) := Wf m+1)T (u) - W5t(uj) for all u £ O, m £ {0, 1, . . . , M - 1} 

M M 

and all M,K £ N. Our numerical scheme which we denote by J-/ 13(H)- 
measurable mappings Y^> M < K : Q -> if at, m £ {0, 1, ... , M}, N,M,K £ N, 
is then given by Y ' ' := Pn(Q and 

Cf* : = ( Y £> M ' K + ^ ■ F(Y^ M ' K ) + B(Y^ K )AW^ K 

~ E ^' W'*) ( 5 M^*) fc) ( 27 ) 
ie^x V J J 

for every m £ {0, 1, . . . , M — 1} and every N, M, K £ N. It can be seen 
that only increments of the Wiener processes W K : [0, T] x Q — >■ [To, -ft' £ 
N, are used in the scheme above and we emphasize that for many SPDEs 
the method (j2"7|) is very easy to simulate and implement (see Sections [Q, @] 
and for several examples) . We now present the convergence result of the 
scheme (ET 



Theorem 1 (Main result). Let Assumptions U§]\ in Section^ be fulfilled. 
Then there is a real number C £ (0, oo) such that 



E 



X™t -Y 

M 



N.M.K 



i 

2 \ 2 



H 



<C\( inf \ \ ^ +( sup fi?\ +M _min(2(7_ ^' 7) ) (28) 

/ioWs /or every m £ {0, 1, . . . , M} and every N, M, K £ N. 

We now explain the result of Theorem [T] more detailed. The root mean 
square difference (E||X™r - Y^' M > K \\ 2 H )^ for m £ {0,1,..., M} and for 
N, M, K £ N of the exact solution of the SPDE (J23]) and of the numeri- 
cal solution (|27|) is estimated in (128]) by a constant times the sum of three 
terms. The first term, i.e. (infj e x\z N \)~^ for iV £ N, arises due to discretiz- 
ing the exact solution spatially, i.e. due to E||JQ — -PatPQHh for iV £ N and 
i £ [0, T]. The second expression, i.e. (sup Jg jr^ K /ij) a for K £ N, occurs due 
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to discretizing the noise spatially, i.e. due to K|| W t — W^Wjj for K E N and 
t E [0, T]. If Uq C U is finite dimensional we choose Jk '■= {j € J\^j 7^ 0} 
for all K E N and obtain (sup je j\j K //j) a = for every K E N in that case. 
The third term, i.e. M~ min ( 2 (T~/ 3 )'T) for M E N, corresponds to the temporal 
discretization error and converges to zero as the number of time steps MgN 
goes to infinity. 



4 Examples 

In this section Theorem [1] is illustrated with various examples. To this end 
let d E {1,2,3} and let H = U = L 2 ((0, l) d , E) be the R-Hilbert space 
of equivalence classes of B((0, l) d )/i3(R)-measurable and Lebesgue square 
integrable functions from (0, l) d to R. As usual we do not distinguish between 
a B((0, l) d )/£>(R)-measurable and Lebesgue square integrable function from 
(0, l) d to R and its equivalence class in H. The scalar product and the norm 
in H and U are given by 



(v, w) H — {v , w) u — j v(x)-w(x)dx 

'(0,l) d 



i 

2 ' 5 



and 

\\ v \\h = \\ v \\u = ( / \v(x)\ z dx 
for all v , w E H = U. Additionally, the notations 



\ v \\c((o,iy,R) '■= SU P \v(x)\e[0,oo] 



and 



v yxj — v(y*) I 

W v \\c r ((oi) d R) := SU P \ v ( x )\+ SU P ~ fi IT — G [0,oo] 

V ' ' ; xG(0,l) d x,y£{0,l) d \\ X ~ y\\R d 

are used troughout this section for all functions v : (0, l) d — > R and all 
r E (0,1]. Here and below we use the Euclidean norms ||x|| R „ := (|xi| 2 + 
. . . + |x„| 2 ) 1/2 for all ) E R" and all n E N. Concerning 

the Wiener process W : [0, T] x Q — > U we assume in this section that the 
eigenfunctions gj E U, j E J , of the covariance operator Q : U — >■ U are 
continuous and satisfy 

SU P lbillc((0,l)d,K) < 00 and (^j ll^llcp((0,l) d ,K)) < °° ( 29 ) 
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for some p G (0,1). We will give some concrete examples for the {gj). & j 
that fulfill (J2HD later. Additionally, we denote by s n (D) G [0, oo), n G N, the 
sequence of characteristic numbers of a compact operator D : H — » H (see, 
e.g., Section 9 in Chapter XI in [9]). Finally, we define the Schatten norm by 



\D\ 



S P (H) - 



X>nOD)M e[o,oo] 



,71=1 



for every compact operator D : H — )■ H and every p G [1, oo) (see also the 
above named reference). 

We now present a prominent example of the linear operator A in 
Assumption 1. Let I = N d and let G H for j Gibe given by 



2 2 sin(ii7rxi) • . . . • sin(^7rx ( i) 



for all x = (xi, . . . , Xd) G (0, l) d and all i = (ii, . . . ,id) G N d . Additionally, 
let k G (0, oo) be a fixed real number and let (Ai) ier be given by 



A; 



K7V 



(< 



«l) 2 + 



for all i = (ii, . . .,id) G N d . Hence, the linear operator A : D(A) C H — >• H 
in Assumption [T] reduces to the Laplacian with Dirichlet boundary conditions 
times the constant k G (0, oo), i.e., 

d 2 \ / «9 2 



w + . . . + 



dx\ 



holds for all t> G -D(A) in this section. Furthermore, let (X/v) 7VgN be given by 
X N = {l,...,N} d for all N G N. 

In order to formulate the drift term in Assumption [2J let /3 = I 



and let / : (0, l) rf x 



— 7- K be a twice continuously differentiable function 

< oo for 



with /(o,i)" \f( x ^°)\ 2(ix < 00 and su Pxe(o,i) d su P?,eR (j^f)( x ,V) 
all n G {1,2}. Then the (in general nonlinear) operator F : Vg — >• H given 
by 

(F(u))(x) = /(a;,i;(a;)) 
for all x G (0, l) d and all v G satisfies Assumption [2] since Vp — Vd C. 

L 5 ((0,l) d ,R) continuously ( see Remark 6.94 in [51] ). 

The formulation of the diffusion term in Assumption [3] is more sub- 
tle. Let b : (0, l) d x 1 -)■ K be a twice continuously differentiable function 
with 



\b(x,0)\<q, 



gn 

dy n ' 



b)(x,y) 



< 



<q (30) 
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and 



J^fcJ (x, y) ■ b(x, y) - \^-b ) {x, z) ■ b(x, z) 



<q\y-z\ (31) 



for all x G (0, l) d , y, z G R, n G {1, 2} and some given g G (0, oo). Then let 
B : Vg — )■ HS(Uq, H) be the (in general nonlinear) operator 



B{v)u ) (x) = 6(x, t>(x)) ■ u(x) 



(32) 



for all x G (0, l) d , t> G Vp and all u G f/o C [/ = H. We now check step by 
step that B : Vg — > HS(U , H) given by (I3"2"j) satisfies Assumption El First 
of all, 5 is well defined. More precisely, we have 

W B ( V )W 2 HS{U ,H) = Yl \\ B ( V )V^9j\\ 2 H = J2^j \\ B (v)9j\\ 2 H 

= y2^\ \b(x,v(x)) ■ gjix)] 2 dx 

- 2^ (J^^ d \Kx,v(x))\ 2 dx ) [ sup \gj(x)\ 2 



and hence 

W B ( V )W 2 HS(Uo,H) 

< 5^M / (\b(x,v(x)) -b(x,0)\ + \b(x,0)\fdx ) \\gj\\ 2 c ( (Q1)dI 

V(0,l) d J V ' ' ' 



< 



3&J 

and finally 



E<^;(7 (K^)i + i) 2 ^)k" 2 

VJ(o,i) d 



C*((0,l) d , 



HS(U ,H) — 1 v\\ V \\h 

<9(II«IIh + 1) (X^'j (j^P INIc((Q,i 
= <lV TT (Q) ( su P IMI C ((o,i) d ,iR)) (NIh + !) < 00 



,i)° 



(33) 
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for every v e Vp which indeed shows that B is well defined. Moreover, B is 
twice continuously Frechet differentiable and we have 

\\ B '( v ) u \\hs(u ,h) = Yl \\( B '( V ) U ) V^9j\\ 2 H = J2^j \\( B '( v ) u )9j\\ 2 H 

jaJ jej 

2 



J2^j 

jej 



(o,i) d 



d 

—b)(x,v(x))-u(x)- gj (x) 



dx 



< y^Vj ( / \u(x) -gj{x)\ 2 dx 
jej \ J (°> 1 ) d 



and hence 



\\ b '( v ) u \\hs(u ,h) < [ J2q 2 N Mil Il*-Ilc((0,i) d , 



\jeJ 



<q\W\\h ( 
\jej / v 



SUP||^|| C ( ( o,i)d )K ) 



= qy/ r Tr{Q) (sup Il^-H 



\j&J 



C((0,l) d ,K) 



u 



H 



for every u, v G Vg which shows 



sup \\B'(v)\\ L{HHS{UoH)) < qy/Tr(Q) (sup ||&|U (01)d R ) 
Additionally, we have 



< oo. 



\\B"(v)(u,w)\\ 2 HS{U(hH) = \\ B "( v )( u i w )\fi r 39jf H 

jej 

dy 2 



jej \ J W d 



b J (x, v(x)) ■ u(x) ■ w(x) ■ gj(x) 



dx 



j&7 V(o,i) d / lv ; 
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and using L 5 ((0, l) d ,R) C L 4 ((0, l) d ,R) continuously shows 
\\B"(v)(u,w)\\ HS{Uo>H) 

\u(x)\ dx 



(0,l) d 



<qVr^Q) (sup|y| 



c((o,i) d 



(0,l) d 



(0,1)" 



tw(x)| cfej ^sup||^-||^ (01)a 



|m(x)| 5 fix 



|w(x)| da; 



(0,1)° 



for every u,v,w G Vg. Therefore, l^Cl 5 ((0, R) continuously shows 

sup \\B"{v)\\ Li2)(y HS{UojH)) < oo 

due to fl29|) and hence, it remains to establish (120]) - (122]) and the symmetry of 
B'{v)B(v) G HS {2 \U .H) for all v G Vp. For the latter one, note that 

[(B'(v)B(v)) (u,u)\ (x) = [B'(v) (B(v)u)u\ (x) 

d 



—b (x, v(x)) ■ b(x, v(x)) ■ u(x) ■ u(x) (34) 
oy J 

holds for all x G (0, l) d , u,u G Uq and all v G Vp which shows that 
B'{v)B{v) G HS( 2 \U , H) is indeed symmetric for all v G Vp. Moreover, 
we have 



\\B\v)B{v) -B'{w)B{w)\\ 2 HS(2)(UQtH) 

= NVk \\B'{v) (B(v)gj) g k - B'{w) (B(w)gj) g k \\ 



< 



/ ( 



d_ 

dy 



b ] (x, v(x)) ■ b(x, v(x)) 



d_ 

dy 



b J {x, w{x)) ■ b(x, w(x)) 



dx \ lls r illa((o,a) d ,R) H5 ,fc llc*((o,i) d ,i;) 



and using (|3T]) yields 

\\B'(v)B(v)-B'(w)B(w)\\ HS(2HUotH) 



< q \\v — w\ 



H ( Yl VjVh j (sup ||^ || 



2 



c((o,i) d 



gTr(Q) (sup H^ll^o,!)",. 



\v — w\ 



H 
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for all v , w E Vp which shows that f[2~Tj) indeed holds. Estimates (TSUI) and 
will be verified in the more concrete examples in Subsections I4.lti4.3l below. 

Concerning the initial value in Assumption |4j let xo : [0, l] d — > R 
be a twice continuously differentiable function with £o|<9(o,i) d = 0. Then the 
•7-o/<6(V^)-measurable mapping £ : Q — > given by £(cu) = xq for all u E Q 
fulfills Assumption H] for all 7 E (0, 1). 

Having constructed examples of Assumptions dHU we now formulate the 
SPDE f[23l) in the setting of this section. More formally, in the setting above 
the SPDE (T53D reduces to 



dXAx) 



KAX t {x) +f(x, X t (x)) dt + b(x, X t (x)) dW t {x) 



(35) 



with X t \Qr Q>1 \d = and X (x) = Xq(x) for t E [0, T] and x E (0, l) d . Moreover, 
we define a family ft : [0, T] x — > R, j E {k E J | fit 7^ 0}, of independent 
standard Brownian motions by 



for all G fi, t G [0,T] and all j E J with /ij 7^ 0. Using this notation, the 
SPDE fl35l) can be written as 



«AXt(x) + /(x,X t (a;)) 



E 



b(x,Xt(x))y/jrjgj{x) dft (36) 



with X, 



1 1 9(0,1)° 



and X (x) = x (x) for t E [0,T] and x E (0, l) d . 



The algorithm (127|) applied to the SPDE (I3"5"j) then reduces to ' M ' 



Pn(xq) and 



N,M,K 
m+1 



Pve A M y 



N,M,K 
m 



T 



+ — . f(. Y N ' 



N,M,K\ 
m 



+ b(.,Y^)-AW m 



M,K 



+ 



2 



. Y 

' m 



■b(;Y T 



N,M,K\ 
m ) 



(37) 

for all m E {0, 1, . . . , M- 1} and all N,M,K E N. Finally, Theorem □ shows 
the existence of a real number C G (0, 00) such that 



E 



(o,i) d 



A T (x) 



< C iV~ 27 + 



SUp /ij 



+ M -min(2( 7 -/3), 7 ) ^ 
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holds for every N, M, K G N. We now illustrate estimate ( 13 8 p in the following 
three more concrete examples. We begin with the introductory example from 
Section □ (see (T7I) and fllgjl). 



4.1 A one-dimensional stochastic reaction diffusion equa- 
tion 

In this subsection let d — 1, T — 1, k — j^, let x : [0, 1] — >■ R be given by 
x (x) = for all x G [0, 1], let f,b : (0, 1) x R ->■ R be given by /(x, y) — 1 — y 
and &(>, y) = jt^ for all a; G (0, 1), y G R, let J" = N, let J K = {1,2,...,K} 
for all Ken, let fij = and let g d = ej for all j G N. The SPDE (JSHD thus 
reduces to 



dXt(x) 



'\2 

X t (x) + 1-X t (x) 



(39) 

with X t (0) = X t (l) = and X (x) = for x G (0,1) and t G [0,1]. The 
SPDE (139|) is nothing else than equation (j!8p in the introduction. In order to 
apply Theorem [T] it remains to verify (12"U1) and (|22|) . Estimate ([2"01 is fulfilled 
for all 5 G (0, 4) here due to Subsection 4.3 in [31]. In order to establish (J22J) 
several preparations are needed. More formally, let (Q, J 7 , P) be a further 
probability space on which a sequence %j : f2 — )• R, 2 G N, of .F/i3(R)- 
measurable independent standard normal random variables is defined. Then 
we define the ^/i3(f/ )-measurable mappings '■ ^ ~ * U by 

X*'V*) :=X)x*(w) (A*)^ei(x) 

i=l 

for all a; G x G (0, 1), K G N and all i? G (0, |). It will be essential to 

estimate E \x K '^{x)\ 2 and E - x K,4 (y)\ 2 for x,y G (0, 1), K G N and 

$ G (0, |) in order to check ([22]) . To this end note that 

K K 



1=1 8=1 

<2(l + ^)(f2i-A (40) 



i=l 
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holds for every x G (0, 1), K G N and every d G (0, \). Moreover, we have 

K 



i=i 



< ^ (A,)- 219 \ ei (x) - ei (y)| 2s (|e,(x)| + \ ei (y)\y 
i=i 

<X>)- W (2vr¥) s 8( 1 -)| a ;- 2/ | 2s 



i=l 

and 

if 



<3(l + /t- 1 ) ^f> (2s_4,9) j |x-2/| 2s 

for every x,y G (0, 1), K G N, i? G (| + \, \) and every s G (0, |). We also 
use the notation 

i 

( f 1 2 f 1 f 1 \v(x) — v(y)\ 2 \ 2 

\\v\\ W r,2 := ( I \v{x)\ 2 dx+ I J — — ^|(i +2r ) dxd V) e[0,oo] 



'0 
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for every £>((0, l))/S(R)-measurable mapping v : (0, 1) — >■ R and every r G 
(0, oo). Then we obtain 



E\\B(v)x 



lvy. 2 ((o,i),K) 
E\b(x,v(x))- X K ''%x)\ 2 dx 

' E \b(x, v(x)) ■ x*'*(x) - b(y, v(y)) ■ X K '*(y) I 




lo Jo \x — y|^ 1+2r ^ 

<2 [ \b(x,v(x))\ 2 E\ x K '*(x)\ 2 dx 



dx dy 



f 1 f 1 \b(x,v(x))\ 2 E\ X K ^(x)-x K ^(y)\\ . 

+ 2 / / 1 — pi dxd y 

Jo Jo \x — y\ 

f 1 f 1 \b(x,v(x)) -b(y,v(y))\ 2 E\ X K >#(y)\ 2 

+ 2 / / 1 — mw) dxd y 

Jo Jo \x — yy ' 

and using f )40|) and fHTj) shows 

<4(l + K - 1 ) |fv«] ||6(-,t;)||^ ((0fl)|R) 



,i=i 

oo 



.1=1 



/o </o \x — y\ 

<4(i+o(E t " (5 ^)ii 6 (-' t; )ii^»((«u)« 



.1=1 



< 

Is — r 



,i=i 
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for every v G H, $ G (f + j, |), s G (r, |), if G N and every r G (0, |' 
Therefore, inequality (23) in Section 4 in [34J gives 



supE\\B{v)x 



KeN 



K,d\\ 2 

I W' 2 ((0,1),R) 



< 



< 



< 



4(i + k- 1 ; 

(s — r) 

4(i + / t - 1 ; 

(s — r) 

12CVg(l + « 
fs — r) 



E^ 2 ^)||&(-^)ll^((o,i), R ) 



(42) 



E 



(2s-4i?) 



i=l 
— 1\ / oo 



(1-r) 



i + IM 



\v\\ Vr ) < oo 



,1=1 



for every # G (| + |, |), s G (r, |), t> G Vr and every r G (0, |). Moreover, 
we have 



-A)~*B{v)Q- 



HS(U ,H) 



(-A)^B(v)Q 
Q^- a) B(v) (-A) 



HS(H) 



-() 



Q 



KIT' 



HS(H) 



B(v) (-Ay 



HS(H) 



B(v) (-Ay 



HS[H,V. i. 



and using inequality (19) in Section 4 in [31] and estimate ( 142]) in this article 
then yields 



;-A)-*B(t;)Q- 



HS(U ,H) 



< 



sup E ||-B(f )x 



Ken 
i\ 2 



K'd\ 



VI) 



<C (a _i) + ( supE||B(z;)x 
12C 2 g(l + K 



K,i9|| 2 

Iw r2 «- 1 . 2 ((0,1),K) 



(a + l- 2a)' 



< 00 



for every 1? G (| + \, \), s G (2a — 1, |), u G V( a -i) an d every a G (|, |). 



Therefore, estimate fl22]) is satisfied for all a G (0, |) and all 7 G (5, §)• This 
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finally shows that Assumptions [T0 are fullfilled for the SPDE ( 15§j) for all 
ae(0,|), /3 = |and all 7 e(±,f). 

Theorem [1] therefore yields the existence of real numbers C r G (0,oo), 



r G (0, |), such that 



•i 4J> 
•1 



E 



o 



X T {x)-Y^ K {x) 



dx 



<C r ( N<r-& + K<r-& + M<r-*> 



(43) 

holds for all N,M,K G N and all arbitrarily small r G (0, |). In order to 
balance the error terms on the right hand side of (1431) we choose N 2 = K 2 = 



M in (j431) and obtain the existence of real numbers C r G (0, oo), r G (0, |), 
such that 



E 



X T (x) - Y, 



N,N ,N , 



iV 2 



X 



dx) <C T - N (r ~^ 



(44) 



holds for all N G N and all arbitrarily small r G (0, |). Estimate (jHJ) 
is nothing else than inequality (|T6|) in the introduction. We also refer to 
Figure [T] in the introduction for numerical results illustrating 



4.2 A one-dimensional stochastic reaction diffusion equa- 
tion with AQ ^ QA 



In Subsection 14 . 1 1 we assumed that the eigenfunctions of the dominating linear 
operator A and of the covariance operator Q of the driving Wiener process 
W : [0, T] x Q — > H of the SPDE (1551) coincide and in particular, we assumed 
in Subsection 14.11 that 

AQv = QAv (45) 

holds for all v G D(A). However, our general setting in Section [2] does not 
need condition (l4"5j) to be fulfilled. To illustrate this fact we consider in this 
subsection an example in which (1431) fails to hold. More formally, in this 
subsection let d = 1, T = 1, k = let x : [0, 1] — > R be given by x (x) = 
for all x G [0, 1], let f,b: (0, 1) X R — >■ M be given by f(x, y) = 1 - y and 
b(x,y) = ^ for all x G (0,1), y G R, let J = {0,1,2,...}, let J K = 
{0, 1, . . . , K} for all K G N, let /i = 0, fij = ± and let g j : (0, 1) -> R be 

given by go(x) = 1, ^(x) = \?2cos(jirx) for all x G (0, 1) and all j G N. The 
SPDE (I35D thus reduces to 



dXAx) 



2 



20 dx 2 



XAx) + 1 - X t (x) 



dt + lfx% dw ^ ^ 
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with X t (0) = X t (l) 
course, the SPDE 



and X (x) = for x 
can also be written as 



G (0,1) and t G [0,1]. Of 



dX t (x) 



20dx* 



XAx) + l-XAx) 



dt 



X t (x) Vg 
+ X t (x) 2 j 1 - 5 



cos(jnx) df3 3 t 



with X t (0) = X t (l) = and X (x) = for x G (0, 1) and t G [0, 1]. Es- 
timate (|2"U|) is here fulfilled for all 5 G (0, |) due to Subsection 4.2 in |34J. 
Moreover, as in Subsection 14.11 it can be shown that inequality (I2"2"j) holds for 
all a G (0, |) and all 7 G (|, 1). This finally shows that Assumptions [TH are 
fullfilled for the SPDE gHD for all a G (0, §), f3 = § and all 7 G (§, 1). 

Theorem [T] therefore yields the existence of real numbers C r G (0, 00), 
r G (0, 1), such that 



E 



X T (x) 



(x) dx) <C r {N^ + K^ + M {r - 1 ^) (47) 



holds for all N,M,K G N and all arbitrarily small r G (0,1). Choosing 
iV 2 = K 2 = M in f T4T|) hence gives the existence of real numbers C r G (0, 00), 
r G (0,2), such that 



E 



dx) <C r - N {r '- 2) 



(4? 



holds for all iV G N and all arbitrarily small r G (0,2). The approximation 
thus converges in the root mean square sense to Xt with order 2— 



Y, 



TVj./V ,7V 
N 2 



as goes to infinity. Since P^{H) C H is iV-dimensional and since time 
steps are used to simulate Y N 2 ' N , 0(N 3 log(N)) computational operations 
and random variables are needed to simulate Y^ N ,N here. Combining the 



computational effort 0(N 3 log(iV)) and the convergence order 2— in 
shows that the algorithm f )3T|) in this article with iV 2 = K 2 = M needs 

3 

about 0(e~~ 2 ) computational operations and random variables to achieve a 
root mean square precision e > 0. 

The linear implicit Euler scheme combined with spectral Galerkin meth- 
ods which we denote by J r /i3(ff)-measurable mappings : Q — > H, 
n G {0, 1, . . . , iV 4 }, N G N, is given by Z$ = and 



7N 

J n+1 



P, 



N 



z 



N 



T 



f(;Z») + b(-,Z»).AW: 



for all n G {0, 1, . . . , iV 4 - 1} and all iV G N here. 



tN 4 ,N 



(49) 
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In Figure H]the root mean square approximation error (E||Xr — Z^ 4 1|^) 2 
of the linear implicit Euler approximation Z^ 4 (see (T4HD) and the root mean 
square approximation error (E||X T — Y^ N ' N \\h)^ °f the approximation 
y^2 N ,N in this article (see (I2"T|) and (13"T|) ) is plotted against the precise num- 
ber of independent standard normal random variables needed to compute 
the corresponding approximation for N G {4, 8, 16, 32}. 




Precise number of used random variables 



Figure 4: SPDE (146]) : Root mean square approximation error (E||Xr — 
Z^4,\\ 2 H )2 of the linear implicit Euler approximation (see fT4"9]) ) and root 
mean square approximation error (E\\X T -Y^ N > N \\ 2 H yioi the approxima- 
tion Y^ N ' N in this article (see fT2T|) and f )3T|) ) against the precise number 
of independent standard normal random variables needed to compute the 
corresponding approximation for N G {4, 8, 16, 32}. 
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4.3 A two-dimensional stochastic heat equation and 
splitting-up approximations 

In comparison to existing algorithms and approximation results in the liter- 
ature the main contribution of this article is to break the computational 
complexity for SPDEs with possibly nonlinear diffusion and drift coeffi- 
cients. However, in order to compare the Milstein type algorithm (1271) in 
this article with previously considered numerical methods we consider a 
linear stochastic heat equation in this subsection. More formally, in this 
subsection let d — 2, T — 1, k — let x : [0, l] 2 — > R be given by 
xo(xi, x 2 ) = 2 sin(7rxi) sin(7nE2) for all X\, x 2 G [0, 1], let /, b : (0, l) 2 x R — > R 
be given by f(xi, X2, y) — and b(xi, X2, y) — y for all x\, x 2 £ (0,1), y G R, 
let J = N 2 , let J K = {1,2,..., K} 2 for all K G N, let fx {juh) = {j x + j 2 )" 4 
and let g(j lt j 2 ) = £(ji,j 2 ) f° r a ^ hi 3i e N. The SPDE (I3"5]) thus reduces to 



dX t (xi,z 2 ) 



50 \dx\ 



d 2 d 2 
+ 



X t (x 1 ,x 2 ) 



dt + X 4 (xi, x 2 ) dW t (x 1 ,x 2 ) 



(50) 

with X t 1 9(o,i) 2 = and X (xi, x 2 ) = 2 sin(7ra;i) sin(7rx 2 ) for x\, x 2 G (0, 1) and 
t G [0, 1]. In view of ([H the SPDE (J50J) can also be written as 



dX t (x 1 ,x 2 ) 



d 2 d 2 
+ 



50 \dxf dxf 



X t (x 1 ,x 2 ) 



dt 



+ E T ■ ( !!'m2 2 sinO'i^i) ^2^2) d/ff lJa) 

with X t 1 9(0,1)2 = and X (xi, X2) = 2 sin(7rxi) sin(7rx 2 ) for x±, x 2 G (0, 1) and 
tG [0,1]. 

Due to Subsection 4.3 in [31] inequality ([2"0~]) holds for all 5 G (0, ~) here. 
In order to verify (I2"2"j) the notation 



L°°((0,1) 2 



:= inf |i2 G [0, c») A({x G (0, l) 2 | v(x) > R}) = o} G [0, 00] 



is used for all £>((0, l) 2 )/i3(R)-measurable mappings v : (0,1) 



in this 



subsection. Then we obtain 

<EK e -^M|e 



L°° ((0,1) 2 ,I 



j|IC((0,l)2,I 



ieN 2 



< 2 E 



-2r 



,ieN 2 




ieN 2 
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for all v G V r and all r G (|, oo). Moreover, we have 



< 



(-A)-*B(t;)Q 
(-^ 



HS(U ,H) 
(i-a) 

H5(i?) 



;-A)-*s(t;)g 



(i-a) 



> ^Hl°°((o,i) 2 



s 



.(H) 



(l-2a) 

(i-a) - 



(l-2oQ 
2 



ieN 2 

and using ( 15Tj) shows 



vjej 



(l-2c) 



J 



(-A)-*S(i;)Q- a 

E ( Ai )~" [ q IMIl-((o,i) 2 ,k) + ?) ( Tr (<2)) 



HS{U ,H) 



< 



(i-a) 



ieN 2 



< 9 (1 + Tr(Q)) (A*)"*) (l + IMU ((0j i )2 , K) ) 

VieN 2 / 

< 2q (1 + Tr(Q)) (A*)"*) f 1 + £ (A,)-* 7 

\i€N 2 / V ieN 2 



i + IM 



< oo 



for every $ G (a, ~), a G (0, A D £ F 7 and every 7 G (|, 1). Inequality fl22|) 
thus holds for all a G (0, |) and all 7 G (|, 1) here. This finally shows that 
Assumptions [Hare fulfilled for the SPDE (JED) for all a G (0, §), /3 = § and 
all 7 G (|,1). 

Theorem [T] therefore yields the existence of real numbers C r G (0, 00), 
r G (0, 1), such that 



1 rl 



E 




Jo 



X T (x 1 ,x 2 ) - Y M ' ' (x 1 ,x 2 ) 



dx\dx'< 



< C r (iV (r ~ 2) + K (r -V + M {r ^) (52) 

holds for all N,M,K G N and all arbitrarily small r G (0, 1). In order to 
balance the error terms on the right hand side of ( )52|) we choose M = N 2 = 
K 2 in (152]) and obtain the existence of real numbers C r G (0, 00), r G (0, 2), 
such that 



E 



fs: 

'0 Jo 



Xt(x\, x 2 ) — Y?fi N ' N (xi, X2) dxidx2) < C r ■ N^ r 2 ' (53) 
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holds for all N £ N and all arbitrarily small r £ (0,2). The approximation 
thus converges in the root mean square sense to X? with order 



N,N 2 ,N 
N 2 



2— as N goes to infinity. The numerical approximations 



N,N 2 ,N . 



n -> h. 



n £ {0,1, 
and 



, iV 2 }, N £ N, (see and (|37j)) are here given by F c 



N,N 2 ,N 







X 



N,N 2 ,N 



n+1 



(54) 



1 + AW, 



N ,N 



+ ^ ( AW, 



v 2 ,v 



2iV 2 ^ 

J'6Jk 



N,N 2 ,N 



for all n £ {0, 1, . . . , iV 2 - 1} and all JV £ N. Since Pjv(#) C H in iV 2 - 
dimensional here and since N 2 time steps are used to simulate Y N i ' N , 
0(N 4 \og(N)) computational operations and random variables are needed to 
simulate Y N 2 ' N . Combining the computational effort 0(N 4 log(iV)) and 
the convergence order 2— in (1531) shows that the algorithm (154|) in this ar- 
ticle needs about 0(e~ 2 ) computational operations and random variables to 
achieve a root mean square precision e > 0. 

The linear implicit Euler scheme combined with spectral Galerkin meth- 
ods which we denote by J-"/i3(if)-measurable mappings : Q — > H, 
n £ {0, 1, ... , N 4 }, N £ N, is given by Z^ = x and 



Pi 



N 



-1 



1 + AW, 



N 4 ,N 



7iV 



(55) 



for all n £ {0, 1, . . . , N 4 - 1} and all JV £ N here. 

Moreover, since the SPDE (l5Uj) is linear here, the splitting-up method 
in [20] can be used in order to solve ( 15"0|) approximatively. The key idea of 
the splitting-up approach is to split the SPDE (150]) into the explicit solvable 
subequations 



dX t (x 1 ,x 2 ) 



1 

50 



d 2 d 2 

+ — 7; ] Xt(xi,X 2 ) 



dx\ 



dx\ 



dt, 



X: 



t\d(o,iy 



= (56) 



and 



dX t (x 1 ,x 2 ) = XAxi, x 2 ) dWAxi, x 2 ) 



(57) 



for t £ [0, 1] and X\, x 2 £ (0, 1). For the solution processes X, X : [0, T] x Q — > 



3 (^-f £ jejW fe) 2 ) 



H of ([56]) and ([57} we obtain X t = e M X and X t = eV" 1 "^^' 
P-a.s. for all t £ [0, 1]. This suggests the splitting-up approximation 



X, 



e At I e 
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for t G [0,1] where X : [0, T] x $7 — >■ H is the solution process of the 
SPDE (|5T3]1 . The resulting splitting-up method which we denote by J r /B(H)- 
measurable mappings Z* : n ->■ #, n G {0, 1, . . . , iV 2 }, JV G N, is then given 
by Zq = xq and 



Z&i = Pn e A ^ \e\ AWn . (58 

for all n 6 {0, 1, ... , A^ 2 — 1} and all JVgN. We remark that I. Gongy and 
N. Krylov considered in [20] temporal splitting-up approximations and we 
added an appropriate spatial approximation here in order to compare the 
splitting-up method with the algorithm (15111 in this article. Using the Talyor 

2 

approximation e x ~ 1 + x + \ for all x G M then yields 

» 1 + - ^5 E + 5 - ip E ft(*> ! ) 

- 1 + AWf> N + \ (AWr*) 2 - 2^ E ^'fe') 2 ( 59 ) 

for all n G {0, 1, . . . , A^ 2 — 1} and all N & N. Using approximation (159]) in 
(158]) finally shows 



7JV 



for all n G {0, 1, ... , A^ 2 — 1} and all A" G N which is nothing else than the 
recursion for Y*> n2 ' n , n G {0, 1, . . . , N 2 }, N G N in ([54]). So, in the case of 
the linear SPDE (150]) an alternative way for deriving the algorithm f[54"|) in this 
article is to apply an appropriate Taylor approximation for the exponential 
function (see (151?]) for details) to the splitting-up method (158"]) . More results 
for the splitting-up method can be found in IS E21 EE US ED 121 EH] and 
the references therein. 

In Figure[5]the root mean square approximation error (E||Xy— Z$± \\ 2 H )^ of 
the linear implicit Euler approximation Z^ 4 (see (155]) ). the root mean square 

approximation error (E\\Xt — Y^ N ' N \\ 2 h)^ of the approximation Y^ N ,N in 
this article (see (154"]) ) and the root mean square approximation error (E||X^ — 
Z^Wh) 1 of the splitting-up approximation Z^ 2 (see (158]) ) is plotted against 
the precise number of independent standard normal random variables needed 
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to compute the corresponding approximation for N G {2,4,8,16,32}: It 
turns out that Z^ 2 4 (32 6 = 1 073 741 824 random variables) in the case of 

the linear implicit Euler scheme ( 15 5 p . that Y^' 32 ' 32 (32 4 = 1 048 576 random 
variables) in the case of the algorithm ( 1541) in this article and that Z^ 2 2 
(32 4 = 1 048 576 random variables) in the case of the splitting-up method (l58j) 
achieve a root mean square precision e = for the SPDE (150|) . The 

Matlab codes for simulating Z|| 4 via d55D, F 3 3 2 2 ' 32 ,32 via (JMD and Z™ 2 via 
(I58p are presented below in Figures El [7] and [8] respectively. The Matlab 
code in Figure [6] requires on an Intel Pentium D a CPU time of about 
29 minutes and 57.06 seconds (1797.06 seconds), the code in Figure [7] 
requires a CPU time of about 1.99 seconds while the code in Figure [S] 
requires a CPU time of about 2.10 seconds to be evaluated on the same 
computer. 

Finally, we conclude that for linear SPDEs the splitting-up method con- 
verges with the same order as the algorithm (12TP in this article. For SPDEs 
with nonlinear diffusion coefficients the splitting-up method can in general 
not be used efficiently anymore since the splitted subequations do in general 
not simplify in comparison to the original considered SPDE. In contrast to 
the splitting-up method the algorithm (127]) works quite well for nonlinear 
equations as demonstrated in Section (TJ Subsection I4.1[ Subsection 14.21 and 
Section 



5 A further numerical example 

Although the setting in Section [2] requires the nonlinear coefficients F and 
B of the SPDE ( I2"3"j) to be globally Lipschitz continuous, we strongly believe 
that our method (127]) produces efficient results for SPDEs with non-globally 
Lipschitz nonlinearities. To substantiate this claim, we apply in this sec- 
tion the algorithm (127]) and the linear implicit Euler scheme to a stochastic 
Burgers equation whose nonlinear drift term F is quadratic and therefore 
not globally Lipschitz continuous anymore. More formally, we consider the 
SPDE 



dXAx) 



-XAx) - XAx) —XAx 
100 dx 2 y ! v ; V dx v ' 



dt + X t (x) dW t (60) 



with X t {0) = X t (l) = and X Q (x) = (sin(Trx) + sin(2vrx)) for x G (0, 1) 
and t G [0, T] on H = L 2 ((0, 1), R) where T = 1 and where W : [0, T] x Q -> 
H is a standard Q- Wiener process with the covariance operator Q : H — )■ H 
given by (Qv) (x) = Xl^Li n ~ 3 sm(mrx) J sin(n7ry) v(y) dy for all x G (0, 1) 
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Precise number of used random variables 



Figure 5: SPDE (|50|) : Root mean square approximation error (E||Xt — 
Z^ 4 ||^)2 of the linear implicit Euler approximation Z^ 4 (see (1531) ). root 
mean square approximation error (E||X^— Y^ N ' N \\%)^ of the approximation 
Ym2 N ,N hi this article (see f l54|) ) and root mean square approximation error 
(E||Xr — Z^ 2 \\ 2 H )^ of the splitting-up approximation Z^ 2 (see ([5S]) ) against 
the precise number of independent standard normal random variables needed 
to compute the corresponding approximation for N G {2, 4, 8, 16, 32}. 
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1 N = 32; M=N"4; Y= zeros(N,N); Y(l,l) = 1; 

2 [nl,n2] = meshgrid ( 1 :N) ; 

s A = -(nl.~2 + n2/2) * pi~2/50; 

4 mu = 1 ./ ( nl . "4 + n2 . A 4) ; 

5 for m=l:M 

e y = dst ( dst ( Y ) ' ) ' * 2; 

7 dW= dst( dst( randn(N,N) .* sqrt(mu/M) )' )' * 2; 
s y = (1 + dW) .* y ; 

9 Y= idst ( idst ( y' )' ) / 2 ./ ( 1 — A/M ) ; 

10 end 

11 surf( nl/(N+l), n2/(N+l), dst ( dst ( Y ) ' )' * 2); 

Figure 6: Matlab code for simulating the linear implicit Euler approxima- 
tion Z% 4 with N = 32 (see (J55J)) for the SPDE (06]). 



1 N=32;M=I\r2;Y = zeros (N,N) ; Y(l ,1) = 1; 

2 [nl,n2] = meshgrid ( 1 :N) ; 

s A = -(nl/2 + n2."2) * pi"2/50; 

4 mu = 1 ./ ( nl . 1 + n2."4); 

s g = zeros (N,N) ; grid = ( 1 :N) / (N+l) ; 

6 for ml=l:N 

7 for rri2=l:N 

8 g = g+4*sin (ml*grid ' * pi ) . " 2 * sin (m2*grid*pi ) . * 2 . *mu(ml,m2) /M; 

9 end 

10 end 

11 for m=l:M 

12 y = dst ( dst ( Y ) ' ) ' * 2; 

13 dW = dst( dst( randn(N,N) .* sqrt(mu/M) )' )' * 2; 
y = (1 +dW+ (dW."2 - g)/2) .* y; 

15 Y = exp( A/M ) .* idst ( idst ( y ' ) ' ) / 2; 

16 end 

17 surf( nl/(N+l), n2/(N+l), dst ( dst ( Y ) ' )' * 2); 

Figure 7: Matlab code for simulating the approximation Y^ N ' N in this 
article with N = 32 (see (J5D) for the SPDE (SSD - 
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1 N = 32; M=N~2; Y= zeros(N,N); Y(l,l) = 1; 

2 [nl,n2] = meshgrid ( 1 :N) ; 

3 A = -(nl. A 2 + n2/2) * pi"2/50; 

4 mu = 1 ./ ( nl . "4 + n2 . " 4) ; 

5 g = zeros(N,N); grid = (1:N)/(N+1); 

6 for ml=l:N 

7 for m2=l:N 

8 g = g+4* sin (ml* grid ' * pi ) . * 2 * sin (m2* gr id * pi ) . * 2 . *mu(ml ,m2) /M; 

9 end 

10 end 

11 for m=l:M 

12 y = dst ( dst ( Y ) ' ) ' * 2; 

13 dW = dst( dst( randn(N,N) .* sqrt (mu/M) )' )' * 2; 

14 y = exp(dW - g/2).*y; 

is Y = exp( A/M ) .* idst ( idst ( y ' ) ' ) / 2; 

i6 end 

it surf( nl/(N+l), n2/(N+l), dst ( dst ( Y ) ' )' * 2); 

Fi gure 8: Matlab code for simulating the splitting- up approximation Z^ 2 
with N = 32 (see (J5BJ) for the SPDE (US}. 



and all v G H. The family gj G H, j G J , of functions with J = N and 
= \/2sin(j7rx) for all x G (0, 1), j G N, thus represents an orthonormal 
basis of eigenf unctions of Q : H — > H. The corresponding eigenvalues (fij)jej 
satisfy \ij = 4 for all j 6 N and we choose Jk = {1,2, ... ,K} for all 
K G N here. Theoretical results for the stochastic Burgers equation including 
existence and uniqueness results can be found in j6l EH EH], for instance. 

The method flZZ} with N 2 = K 2 = M applied to the SPDE (JSUj) is given 
by If'^ = P N (X ) and 

N \ 

+ y^' n2 ' n ■ (AWf> N + \ (AWf> N ) 2 - ^ J> fe') 2 ) ( 61 ) 

3=1 / 

for all n £ {0,1, ... , N 2 — 1} and all iV G N here. 

The linear implicit Euler scheme applied to the SPDE (16 Op which we 
denote by J 7 / B(H)-mea.smab\e mappings : ft — > H, n G {0, 1, . . . , N 3 }, 
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iV G N, is given by Z$ = P N (X ) and 

d = ^ ^ ^ -^- Z n- (§-X Z ") + Z " • AW " 3,N ) (62) 

for all ne {0,1,..., N 3 — 1} and all iV G N here. 

Recently, it has been shown in [28] that Euler's method fails to converge in 
the root mean square sense to the exact solution of a SODE with a superlin- 
early growing drift coefficient of the form (160]) (see Theorem 2 in [28] for the 
precise statement of this result) although pathwise convergence often holds 
(see, e.g., HUES] for SODEs and HB3 SS] for SPDEs). That is the reason 
why we consider the pathwise instead of the root mean square approximation 
error in this section. More precisely, in Figure [9] the pathwise approximation 
error ||Xr(o;) — Z^ 3 (u)\\ 2 H of the linear implicit Euler approximation Z^ 3 
(see fl62l ) and the pathwise approximation error || Xr(u;) — Y^ N ' n (oj)\\h of 
the approximation Y^f 1 ' N in this article (see (16T]) ) is plotted against the 
precise number of independent standard normal random variables needed to 
compute the corresponding approximation for iV G {2,4,8,16,32,64,128} 
and one random u G Q. 



6 Proof of Theorem H 

The notation 



\Z\\ LP{ n;E)--=(^\\Z\\ P E y G[0,oo] 



is used throughout this section for an M-Banach space (E, a J '/8(E)- 

measurable mapping Z : Q — > E and a real number p G [1, oo). We also use 
the following simple lemma (see, e.g., Lemma 1 in |34J and also Theorem 37.5 
in [56]). 

Lemma 1. Let Assumptions in Section^ be fulfilled. Then we have 

IK-^) r ^IL(H) < 1 and U-^r (*" - i)\\hb) ^ 1 

for every t G (0, oo) and every r G [0, 1]. 

In the following Theorem [T]is established. First of all, note that the exact 
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Precise number of used random variables 



Figure 9: SPDE fl60|) : Pathwise approximation error — Z^ z (u))\\ 2 H of 

the linear implicit Euler approximation Z^ 3 (see (16"2"j) ) and pathwise approx- 
imation error \\X T (u) - Y^ N (oo)\\ H of the approximation Y^ N in this 
article (see (16 ip ) against the precise number of independent standard normal 
random variables needed to compute the corresponding approximation for 
N G {2, 4, 8, 16, 32, 64, 128} and one random u G fl 
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solution of the SPDE f[2"31 satisfies 



X mh = e Amh i+ / e A{mh - s) F{X s )ds+ / e A{mh - s) B(X S ) dW s 
Jo Jo 

m-l n(l+l)h m ~ l r(l+l)h 

e A(mh - s) F(X s )ds + J2 / e A{mh - s) B{X s )dW s 

;_n J lh in J lh 



^Amh 



1=0 " " l 1=0 



(63) 



P-a.s. for every m G {0, 1, . . . , M} and every M G N. Here and below h is 
the time stepsize h = Hm — T7 with Mel In particular, (163]) shows 

z=o / 

/m-l /.(J+l)/, \ 

+ PN [Y,J lh e Ai - mh -^B{X s )dW s \ (64) 



P-a.s. for every m G {0, 1, . . . , M} and every N, M G N. In order to estimate 
the difference of the exact solution (163]) and the numerical solution (157)1 we 
rewrite the numerical method ( l2Tj) in some sense. More precisely, the identity 



(Y£' M ' K ) (B(Y^ K )AW^ K y Wl 



M,K 
m 



1 E*#( Y ^*)( B W' M ' K )9 i )9 j 



2M 



(m+l)T 

M 



B , (y wj ^ jT 5 (y W) ^ ^ (65) 



mT 

A/ x " M 



P-a.s. holds for all m G {0, 1, . . . , M - 1} and all N,M,K G N. The proof 
of (165|) can be found in Subsection 16.71 Using (165|) shows that the numerical 
solution (EZD fulfills 



(m + l)T 

.V..W.M dW K 



(m+l)T \ 

+ / " 5' (V™) ( f ^ ( Y™) dW?) dW* (66) 
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P-a.s. for every m G {0, 1, . . . , M — 1} and every N, M, K G N. Therefore, 
the numerical solution (|27p satisfies 



/m-l 

\ z=o ™ V 

(m— 1 \ 
W e^" 1 -^^^)^ (67) 

tm— 1 ^(i+l)^ / r s \ ^ 

£ jf e A (— 0*5' (>f (J B (Y t N ' M ' K ) dW«) dWf 

P-a.s. for every m G {0, 1, ... , M} and every N, M, K G N. In order to ap- 
praise E ||X m/l - Y^ M ' K \\ 2 H for m G {0, 1, ... , M} and N,M,K G N, we de- 
fine the J-"/i3(if)-measurable mappings Z^' M,K : Q — > H, m G {0, 1, . . . , M}, 
N,M,K E N, by 



V / ^"^(X^cW. 
z=o 

/m-l .(f+i)fc / „ s \ 

+ Pn[T / e A ^ h B'(X lh ) / fl(X, A ) <W* dW t 

\ l=Q Jlh \Jlh J 



(6* 



P-a.s. for every m G {0, 1, . . . , M} and every N, M, K G N. The inequality 

(ai + . . . + a n f < n ((ai) 2 + . . . + (a n ) 2 ) (69) 
for every a\, . . . , a n G R and every nGN then shows 

E ||-X" m h — ^m^'^^ll^ < 3 • E ||X m/i — P A r(X m / l )|| 2 ^ 

+ 3 ■ E \\P N (X mh ) - Z»> M ' K \\ 2 H + 3 • E ||Z^ - Y»' M > K \\ 2 H (70) 

for every m G {0,1,..., M} and every N,M,K G N. In order to esti- 
mate the expressions E \\X mh - P N (X mh )\\ 2 H , E \\P N (X mh ) - Z% M ' K ||^ and 

E \\Z% M > K - Y^ M ' K \\ 2 H for m G {0, 1, ... , M} and N,M,K G N, the real 
number i? G (0, oo) satisfying 

E ||fl(X t )||L(i*,,v.) < ^ \\F'(v)\\ L(H) < R, \\F"(v)\\ Li2)(v ^ H) < R, 
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E \\ F { X t)\\ 2 H < R , \\ B '( v )\\l(H,HS(U ,H)) ^ R > \\ B "( V )WlW(V ,HS(U o ,H)) ^ R i 

E\\(-AyX t \\ 2 H = E\\X t \\ 2 Vi < R, E\\X t2 -X tl f v , < R \h - t^^'^ , 



c + 



1 



1 



1 



+ T+ A 



-H 



\L(H) 



< R 



(1-7) (1-20) (1-25) 

for every v G Vg and every t, ti , t 2 G [0,T] is used throughout this proof. 
Due to Assumptions dH in Section [2] and Proposition [TJ such a real number 
indeed exists. For the spatial discretization error E \\X mh — i"jv(-^m/i)||jj we 
then obtain 



EIIX 



mh 



Pn{X, 



mh) \\h 



E\\(I-P N )X,, 



N) s*-mh\\H 



E 



-Ay 1 (i-p N )(-Ayx,, 



mh 



H 



(71) 



<\U-A)-^I-P N )\ 



L{H) 



EIIX 



mh 1 1 y 



< R (r N ) 



for every m G {0, 1, . . . , M} and every N,M G N where here and below the 
real numbers (r/v)jveN c ^ are gi ven by 



r N :-- 



-A)-' (I-P N )\ 



L(H) 



inf Aj 



for every iV G N. The rest of this proof is then divided into six parts. In the 
first part (see Subsection 16.11) we establish 



E 



TO — 1 

E 

1=0 



(l+l)h 



D A(mh—s) 



F(X S ) - e 



_ A(m—l)h 



Ih 



F{X lh )) ds 



< 



3QR 8 



jVfmin(4( 7 -/3),2 7 ) 



H 



for every m G {0, 1, . . . , M} and every M G N. We show 

2 



E 



m-l p(l+l)h 

/ e A{mh - s) B(X s )d(W s -W s K ) 



< 4R 5 



H 



SUp fJLj 

jeJ\J K 



2a 



(72) 



(73) 



for every m G {0, 1, . . . , M} and every M, K G N in the second part (see 
Subsection 16.21) and 



E 



m— 1 

E 

1=0 



(l+l)h 



D A(mh-s) A(m— l)h 



Ih 



) B{X S ) dW : 



K 



< 



3R 4 



H 



M( 1+25 ) 



(74) 
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for every m G {0,1,..., M} and every M,K G N in the third part (see 
Subsection 16.31) . The fourth part (see Subsection I6.4p gives 



E 



J2 / e A ^ h / B"(X lh + r(X s - X lh )) (X s - X lh , X s - X lh ) 
1=0 Jlh Jo 



■ (l-r)drdW* 
and in the fifth part (see Subsection I6.5P we obtain 

m ~ l n(l+l)h 



< 



R 6 



j^min(4( 7 -/3),2) 



(75) 



H 



E 



J2 / e A ^ h B'(X lh ) (X s -X lh - B{X lh ) dW« dW s 

,_ n Jlh V Jlh J 

+ 20R n ( sup n^) (76) 



j\//min(4(7-/3),2 7 ) 



for every m G {0, 1, . . . , M} and every M, if eN. The inequalities (JZ2J)- dZS]) 
are used below to estimate E [^^(X^) — Z^ M ' K \\ 2 H for m G {0, 1, ... , M} 
and N, M, K G N in (T70|) . In the sixth part (see Subsection 16.61) we estimate 



E\\Z. 



N,M,K _ y N , M ,K\\ 2 



9R 



'm—l 



N,M,K 



1=0 



(77) 



for every m G {0, 1, ... , M} and every N, M, K G N by using the global 
Lipschitz continuity of the coefficients F : Vp — > H (see Assumption [2]) and 
B : V p -> HS(U ,H) (see Assumption EJ. Combining (170]) , (171]) and (177j) 
then yields 



EX 



1". 



Af,M,_ft:|| 2 



<3i?(r w ) 5 



+ 3-E P^X^)-^' 



N,M,K 1 1 2 . 27-R 4 



'm— 1 



N,M,K 



1=0 



for every m G {0, 1, ... , M} and every X, M, A" G N. Hence, fH]), fl68]) and 
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show 



ME 



. 1=0 



w\\y vN,M,K 
< 9-E 

+ 9-E 
+ 9-E 



i2 
111 



m-l r (l+l)h 



I Hi 



{e Mmh - s) F(X s ) - e A{m - l)h F(X lh )) ds 



H 



ME 



™-l r (l+l)h 



, 1=0 
' m— 1 



a A(mh— s) 



ME 



«=0 



c(i+l)ft 



S(X s )d(W s - wf) 



H 



A(mh—s) 



B(X S ) - e A ^ h B(X lh )) dW t 



K 



m— 1 

E 



a A(m- 



' 1 >B\X 



lh, 



Ih 



B(X lh )dW«)dW. 



lit 



rK 



H 



27R 



f rn—l 



+ 3R(r N ) 2 + —[J2^\\Xi h -Y l 



N,M.K 



. 1=0 



and using ||-Pjv(iOll.ff < \\v\\ H for all v & H additionally gives 



E\\X mh -Y, 



N,M,K 
m 



|2 
111 



< 9 • E 
+ 9-E 
+ 9-E 



m-l r (l+l)h 



lh 



E 



m— 1 

E 

1=0 



m -! r (l+l)h 



A(mh—s) 



a A(m-l)h 



F(X lh )) ds 



(l+l)h 



A(mh—s) 



lh 



B(X s )d(W s -W s K ) 



H 



E 



(e A{mh - s) B(X s ) - e A{m - l)h B(X lh )) dW : 



K 



m -! ?{l+l)h / rs \ 

J2 / e A ^ h B'(X lh ) / B(X lh ) dW* dWf 

,_ n Jlh \Jlh J 



+ 3R (r N y 



1=0 

27R 4 



H 



f m— 1 



M 



E E \\ X ^- Y i 



N,M,K 



. 1=0 



for every m G {0, 1, . . . , M} and every N, M, K G N. Therefore, f[?2]) and 
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flTBJ yield 



E 1 1 X mh — Y r 



N,M,K\\ 2 



< 



324R 8 



+ 18 ■ E 
+ 18 • E 



m OH — j^min(4(7-/3),27) 

Ih 



+ 36 i? 5 I sup /ij J 

\jeJ\J K J 



2a 



E 

1=0 

m—l 

E 



D A(mh— s) A(m—l)h 



K 



1=0 

m -i "(l+l)h 



H 



A{m-l)h 



{B{X s )-B{X lh ))dW t 



K 



11, 



m—l 

E 



1=0 



(l+l)h 



e A{m - l)h B'(X, 



u, ) 



a. 



B{X lh )dW*\dW? 



111 



H 



27R A 



'm—l 



+ 3R(rN) 2 + ^-[J2 E \\ X ^- Y i 



N,M,K 



1=0 



and f[T4"j) shows 



E \\X mh — Y„ 



N,M,K\\ 2 



< 



32AR 8 



771 OH — ]^min(4(7-/3),27) 



3QR° 1 sup /ij J + 
\jeJ\J K J 



2a 



MO+25)" 



+ 18-E 



E ' 



1=0 



ih 



m—l 



-E 



(l+l)h 



Hi 



e A(m ^ h (B(X s )-B(X lh ))dW s K 

l)h B\X lh ) ^B{X lh )dW^j W\ 



A(m- 



H 



'm—l 



rN.M.K 



1=0 



for every m e {0, 1, . . . , M} and every N, M, K £ M. The fact 

BpQ - S(Jf, fc ) = B'(X lh )(X s - X lh ) 

+ / B"(X lh + r(X s - X lh )){X s - X lh , X s - X lh )(l - r) dr 
Jo 
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for every s G [Ih, (l + l)h], I G {0, 1, . . . , M - 1} and every MeN then yields 

?8 / \ 2a 



E||X m ,-YjW|| 2 „< 



324i? 8 



+ 36-E 
+ 36-E 



ff - j^min(4( 7 -/3),2 7 ) 



54i? 4 



]Vf(l+25) 



T / e A (— V h B\X lh ) (X s -X lh - B{X lh ) dW*\ dW t 

l=0 Jlh V Jlh J 

V / e A ^ h / B"(X lh + r(X s - X lh )) (X s - X lh , 

X s -X lh ) (l-r)drdW s K 



H 



H 



'm—1 



M,K 



1=0 



for every m G {0, 1, . . . , M — 1} and every N, M, K G N. Therefore, ([75 
and (1751) give 



E II y h - 



< 



324# 8 



+ 756i? n ( sup fij ) + 



// - j^min(4( 7 -/3),2 7 ) 

720i? 13 

+ 



+ 



jeJ\J K 
36R 6 



2a 



54fl 4 

M( 1+2<5 ) 



j^min(4(7-/3),2 7 ) j^min(4( 7 -/3),2) 



+ 3i?M 2 + ^ 



^m— 1 



M,K 



, 1=0 



and hence 



E X mh -Y 



N,M,K\ 
m 



2 H < (32AR 8 + 54# 4 + 720i? 13 + 36.R 6 ) 



j£J\Jk 



2a 



27i? 4 



/- m— 1 



+ 756i? n ( sup fij ) + 3i? (r.v) 2 + — ( ^ E \\X lh - Y, 



^-min(4( 7 -/3),2 7 ) 
2 



rN,M,K 



1=0 



H 



for every m G {0, 1, . . . , M} and every iV, M, if eN. Gronwall's lemma thus 
shows 



^\\X mh -Y^ K \\ 2 H <e^ 



( (324i? 8 + 54i? 4 + 720i? 13 + 36i? 6 ) 



j^min(4( 7 -/3),2 7 ) 



+ 756i? ii ( sup fij) +3R(r N y 

jeJ\J K 



2a 
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and hence 



E \\X mh -V 



N,M,K\\ 2 



m 1 1 // 



< 1134i? 13 e 27ij4 



2a 



(r N y+[ sup fij] + M- min(4( ^- /3) - 27) 
jeJ\J K 



for every m G {0, 1, . . . , M} and every N, M, if el. Finally, we obtain 



- m ||// 



< 34R 7 e 14Ri 



inf X-\ ( sup fMj] -\ jr n,uu2r - 
iei\i N J \jej\j K 



< e 



20R 4 



inf Xi) +[ sup ft) + M" min(2(7 - /3) ' 7) 

i&\X N J \j£j\J K 



for every m G {0, 1, . . . , M} and every N, M, K G M. 

6.1 Temporal discretization error: Proof of ( 172| ) 

First of all, we have 



/ (e A(m?i - s) F(X s ) - e A(m - /)/l F(X^)) ds 



< 



m-l f (i+i)h 

£ / IK- 



1=0 



A(mh-s) A(m-l)h 



+ 



Ih 

m ~! r (l+l)h 



)F(X S 



\L 2 (n-,H) 



L 2 (Q;H) 

ds 



2=0 



A(m-l)h 



{F{X s )-F{X lh ))ds 



ih 



L 2 (n-H) 



for every m G {0, 1, . . . , M} and every M G N. Using 

F(X S ) - F(X lh ) = F\X lh ) (X s - X lh ) 

+ / F"(X lh + r(X s - X lh )) (X s - X lh , X s - X lh ) (1 - r) dr 
Jo 
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for all s E [Ih, (I + l)h], I e {0, 1, . . . , M - 1} and all M e N then shows 



£ / (e A ( mfe - s )F(X s ) - e A ( m - ; ) /l F(X^)) ds 



Z=0 

m ~ l r (l+l)h 



Z=0 



\^A(mh— s) gA(m— Z)/i| 



m— 1 



L 2 {Q;H) 
L(H) \\ F ( X s)\\v(Cl;H) ds 



p(l+l)h 

£ / e A ( m -^F'(X^) (X s - ds 
z=o ^ 

"<■-! r(l+l)h nl 

+ J2 / ||F // (^ + r(X s -X^))(X s -X^X s -X„ i )|| LW) 



L 2 (Q;H) 



(1 — r) <ir (is 



and hence 



m— 1 

E 



(H-i)fc 



A(mh— s) 



Ih 



F{X S ) - e A ^ h F(X lh )) ds 



L 2 (fl;H) 



~(l+l)h 



< R\2h + 



\ A(mh—s) _ A(m—l)h\ 



lib A n 

>•*-*■ p(l+l)h 

/ e A ( m -^F'(X^)(X s -X^)rfs 
z=o ^ 

<"-l r(l + l)h pi 



■ ■ 2 



L 2 (^;R) 



L 2 (Q;H) 

(1 — r) ds 



for every me{0,l,...,M} and every MgI Therefore, we obtain 



/ (e A(mft - s) F(X s ) - e A(m ~ l)h F(X lh )) ds 



L 2 (Q;H) 



m-2 



<R[2h + J2 



(l+l)h 



+ 



1=0 

m-1 ~(l+l)h 



Ih 



\\ Ae A(mh- 



s) I 



\L(H) 



| A- 1 ( e ^-^)_/ 



L(H) 



ds 



£ / e A (-^F'(^,) (X. - X Ifc ) ds 



L 2 (n-,H) 



R 



f m— 1 



+ T E 



//i 



E||X S -X, 



|4 \ 2 



(is 
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and Lemma [T] gives 

m ~ 1 r (l+l)h 



/ (e A{mh - s) F(X s ) - e A{m ~ l)h F(X lh )) ds 
1=0 J lh 



L 2 (fl;H) 



'Hi 



<Rl2h + J2 

V 1=0 
m -i r (l+l)h 

+ E 

1=0 

p(l+l)h 

+ E/ 

1=0 Jlh 
™—l r (l+l)h 

E 



m -2 p(l+l)h 



'Ih 



(s - lh) 

(mh — s) 



ds 



a A(m- 



l)h F\X, 



111 I 



+ 



e A(m-l)h pi 



Ih ) 



1=0 



Ilk 



I e A{s - u) F(X u )du) ds 

'lh J L 2 (Q;H) 

e A{s - u) B(X u ) dwA ds 



'lh 



L 2 (Q.;H) 



+ ?(E 



' m - 1 /•(i+l)h 



. 1=0 



llh 



R( S -ihr m( ^~^y ds 



for every m G {0, 1, . . . , M} and every M e N. This shows 

r(l+l)h 

hh 

^ f^ h (s-lh) 



m - 1 r(l+l)h 

/ (e A(mh - s) F(X s ) - e A(m ~ i)?i F(X^)) ds 
i=o J lh 



L 2 (Sl;H) 



(m-2 Q 
1=0 Jlh 



(m — I — l)h 



ds 



m ~ l f(l+l)h 



+ E 



1=0 



m -! r (l+l)h 

+ E 



1=0 
m—l 



lh 



\F'(X lh )((e 



F\X lh ) ( f 

\Jlh 



A(s-lh) 



I) x lh ) 



ds 



e A{s - u) F(X u ) du 



ds 



L 2 (fl;H) 



+ { E E 



1=0 



r {l + l)h / ns \ 

/ e A ^ h F\X lh ) / e A ^B(X u ) dW u ds 

Jlh \Jlh / 



2^2 



H 



d2 A™" 1 

+ ( Yl /i (i+min(2(7 ' /3) ' i)) 



. z=o 
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and 



m ~ 1 r(l+l)h 

/ (e A{mh - s) F(X s ) - e A{m ~ l)h F(X lh )) ds 
1=0 ^ lh 



L 2 (Q;H) 



m-2 



<R[2h+J2 



h 



+ i jR 2 T/l mi„(2( 7 - / 3),l) 

^2(m- 1-1) ) 2 

;=o v 'J 



(m-\ r a+i)h 

z=o 



e/ / 

E/ e 

l_n Jlh 



e Ms - lh) -l)X lh \\ L2{n . H) ds 
e Ms ' u ^F(X u ) du ds 



L 2 (U;H) 



F\X lh ) 



,A(s-u) 



B(X u )dW u ) 



Ih 



ds 



H 



for every to e {0, 1, ... , M} and every MgI Hence, we have 



m -! r(l+l)h 

/ (e A(m ^ s) F(X s ) - e A(m - z)/l F(X^)) ds 
i=o ^ lh 



L 2 (Q;H) 



< 



5 min(2( 7 -/3),l) 



-Ar (e A{s - lh) -l)\\ L{H) \\(-AT X lh \\ L2 ^ H) 
\\ F ( X u)\\ L 2 (n . H) duds 

+ RVh\J2 E / e A(s - u) 5(X u )ciH4 



r {l+l)h 

E 

E/ , 

( m-1 f (l+l)h 



ds 
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and Lemma [T] shows 



/ {e A( - mh ~ s) F(X s ) - e A( - m ~ l)h F(X lh )) ds 
<R^2h+^(l + ln(M))^ + i^^-m^CT-^),!) 



L 2 (Q;H) 



'm— 1 



(Z+l)/i 



V 1=0 

/ E||e A(s - M) 5(X n 



HS(Uo,H) 



du ds 



for every m G {0, 1, . . . , M} and every M G N. Therefore, we obtain 



m—l 



(l+l)h 



i=0 



A(mh— s) 



F(X S ) - e A ^ h F(X lh )) ds 



L 2 (Q;H) 



< Rh f^ + l i n (M)) + -R 4 M- min ( 2 ^-/ ? ). 1 ) + R 2 Mh^ + -R 3 h 
V 2 2 y 2 2 

f m - 1 p(l+l)h rs 
/ / ®\\B(X U )\\] 



\HS(U ,H) 



du ds 



and hence 



/ (e A(m/l " s) F(X s ) - e A{m - l)h F{X lh )) ds 
i=o J lh 

C R 2 M~ l (- + - ln(M) > ) + i jR 4 M -min(2( 7 -«,l) + + I^M" 1 

\ 2 2 / 2 2 



L 2 (U;H) 



{ m -! [>(l+l)h ps 
E / / M) 
i=0 </z/i ,/z/l 



E 



HS(U ,H) 



du ds 



for every m G {0, 1, . . . , M} and every M G N. This yields 



L 2 (Q;H) 



m ~ 1 r(l+l)h 

/ (e A{mh - s) F(X s ) - e Mm ~ l)h F(X lh )) ds 
i=o J lh 

< %RL 2 M- X (1 + ln(M)) + 1^^-^(2(7-^),!) + j£ M ~i + -tfM- 1 

+ R 2 Vh(^MhA 
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for every m G {0, 1, . . . , M} and every M G N. The estimate 



1 + ln(x) = 14 / - ds < 1 + 



ds = 1 + 



8=1 



^ (a; 1 * — 1) x r (1 — r) ^ 

fy iy 



for every r G (0, 1] and every x G [1, oo) then shows 



m-l n(l+l)h 

/ (e A{mh ~ s) F(X s ) - e A{m - l)h F(X lh )) ds 



2=0 

<5 tf M(1 " 7) 



L 2 (Q;H) 



+ 



i? 4 



2 M(l- 7 ) 2M min ( 2 (T-«' 1 ) M 7 2M 



5i? 4 
< — + 



R 4 



R 4 R 4 R 4 
and finally 



m-l p(l+l)h 

/ (e A{mh - s) F{X s )-e A{m ~ l)h F{X lh ))ds 



L 2 (U;H) 



R 4 



< 



6R 4 



5 1 1 

2 + 2 + 1+ 2 +1 / M min ( 2 (T-^).7) ~ Mmin(2( 7 -/3), 7 ) 



for every m G {0, 1, . . . , M } and every M G N. 



6.2 Noise discretization error: Proof of ( 1751 ) 

In this subsection we have 

2 



E 



J e A ^B(X u )d(W u -WK) 



H 



E 



f e A{t - u) B(X u )g 3 d(g^W u ) l 



3&J\Jk 



H 



to ( [ E\\e A ^B(X u ) gj \\ 2 H du) 

<T\ T„ \J S J 
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and hence 



e Mt-v) B ( Xu )d(W u -W?) 
jeJ\J K KJs 



jeJ\J K 



for every s, £ G [0, T] with s < t and every ifeN. This shows 



E 



j\ A ^B(X u )d(W u -W«) 



H 



<( 



SUp flj 

\jej\j K 



2a 



£ ^ ( rE||^- u )s(x u )g- a ^|^d« 

. ^ 7 \./s 



, jeJ\J K 



( 



< I SUp fj,j 

\jeJ\J, 



7 Vie^ Js 



and 



E 



e A ^B(X u )d(W u -W«) 

2a / «t 



ii 



lfrS(t/o,#) 



HS(U Q ,H) 



du 



<( sup ^ (/ E||e A (*- u ) J B(X u )g^ 
\jeJ\J K J \Js 

<( sup ^ f /* (t-«)- M E (-A)-^(x u )g- a 

\jeJ\J K J \Js 
for every s, t G [0, T] with s < i and every if eN. Therefore, we obtain 

nt 2 

E / e^- u )5(X tt ) d(W u - W^) 

<c 2 ( sup /i^ (/ (t-^-^E^fl + HXj 
\jeJ\J K / \Js L v 



V- 



<2c 2 (^ sup /i^ Qf (t-?/) -2 * (l + E||X M ||^) du 
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and 



E 



e A ^B(X u )d(W u -W«) 

2a f r t 



H 



< AR 6 ( SUp fij 

J<eJ\Jk 



(t — u) 2 ^ du 



AR [ sup fij 

jeJ\J K 



2a / f (t-s) 



u 2 Uu 



for every s, t £ [0, T] with s < t and every K £ N. Hence, we have 



E 



I e A ^B(X u )d(W u -W^) 

J s 

< AR 3 1 sup 



2a r u (l-2tf) 



u=(t— s) 



u=0 



< AR A ( sup /ij J it - s) 



(l-2i?) 



and finally 



E 



f e A ^B{X u )d(W u -W«) 

J s 



2 / X 2a 

5 i 



H 



< AR sup (7, 



for every s, t £ [0, T] with s < t and every if 6N. In particular, we obtain 

2 



E 



/ e A{mh - s) B(X s )d(W s -W s K ) 



E 



mil 



i A{ - mh - s ^B{X s ) d(W s - W S K ) 



2 / \ 2c 

H 



< AR b sup fij 



jeJ\J K 



for every m £ {0, 1, . . . , M} and every M, K £ N which shows ([73 
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6.3 Temporal discretization error: Proof of ( 1731 ) 

Here we have 



E 



TO— 1 



(Z+l)/i 



£ . ( 



,A(mh—s) _ A(m—l)h 



1=0 " lh 
m ~l '(l+l)h 



)B(X S ) dW t 



K 



H 



<J2 E II (e A(m/i " s) - e Mm - l)h ) B(X S ) | 



HS(U ,H) 



ds 



1=0 

TO-l /.(/+!)/! 



<V/ (_A)- 5 ( e ^ mh - s ) -e A(m " E (— 

, n Jlh L ( H ) 



1=0 

and hence 



HS(U ,H) 



E 



TO— 1 



(l+l)h 



A(mh—s) A(m—l)h 



) B(X S ) dW t 



K 



t™- 1 Ml+l)h 
J2 / (- A )- S ( e Mmh-s) _ e A(m-l)hj 

,-n Jlh 



< R 



1=0 

rah 



H 

2 

L(H) 



ds 



W_Ayt ^ p A(mh-s) _ p Ah\ 



(m-l)h 
m ~ 2 »(i+l)ft 



ds 



E / ik-a)- 1 ^-^-/) 



|2 

li(H) 



(1-5) A(mh-s) 



ds 



L(H) 



for every m G {0,1,. . . , M} and every M, KeN. Therefore, we obtain 

2 



E 



m— 1 



1=0 



(l+l)h 



£ . 



A(mh— s) A(m— l)h 



< R 



lh 
mh 

(m-l)h 



)B{X s )dW ! 



H 



(m-2 
£/ 



L(ff) 



_ A ^(lS) e A(mh-s) 



L(H) 



ds 
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and 



E 



m ~ 1 r(l+l)h 

V / [e^mh-s) _ e A[m-l)h^ B ( X j dw l 

,_n Jlh 



< R 



Ih 

mh 
(m-l)h 



(s - (m - l)hf 5 ds + Rh 



H 

2 r (i+i)h 



1=0 Jlh 



(mh — s) 2 ^ 5 ^ ds 



'm— 2 



< Rh^ + Rh 3 lJ2( m - 1 - h2{6 ~ 1} 

V 1=0 

for every m G {0, 1, . . . , M} and every M, K G N. This shows 

2 



E 



m -! r (l+l)h 

£/ ( 



3 A(mfe-s) _ A(m-l)h 



H 



f m—1 



< RhP +7S > + Rh (1+2 ^ \£ l2[S ' 1] ) < R h(1+25) [ 2 + ^/ 2{ ^ 1} 



and finally 



E 



m— 1 

E 



(■ 



3 A(m/i-s) _ A(m—l)h 



K 



H 



< Rh^ \2 + 

< 3R 2 h^ < 



■ s (2«-l) 

.(2^1) 
3i? 4 



s=l 



(1 - 25) 



M(l+25) 

for every m G {0, 1, . . . , M} and every M, K G N. 
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6.4 Temporal discretization error: Proof of 



We have 

m— 1 



E 



m-l p(l+l)h pi 

i=0 Jlh Jo 



(1 -r) rfrdW ; 



A' 



H 



m ~ l p(l+l)hpl 

<J2 / E\\B"(X lh + r(X s -X lh ))(X s -X lh ,X s -X t 

l=Q Jlh JO 

for every m G {0,1,. . . , M} and every M, K G N. Therefore, we obtain 

m-l p{l+l)h pi 

J2 / e A ( m "^ / fl"(X„, + r (X s - (X, - X lh , X s - X lh ) 

,_n Jlh JO 



■ (l-r)drdW ! 



K 



H 



m-l r Q+l)h 



lh 



R Xc — X, 



iii 



(m-l „(i 
l=0 Jlh 



-(l+l)h 



E\\X X -X lh \l, ds 



and hence 

m— 1 



E 



m-l p(l+l)h pi 

V / e A ^ h / S"(X Ih + r (X s - (X s - X lh , X s - X lh ) 
i=o Jlh J ° 



• (1 -r) drdW ! 



A 



H 



1=0 



'm—1 



'm-l f (l+i)h 



«=0 



i? 6 



< jR 3 M ^(l + min(4( 7 -«,2)) = jR 3 T ^min(4( 7 -«,2) < _ 
- - y^fmm(4( 7 -/3),2) 

for every m G {0, 1, . . . , M} and every M, G N. 

6.5 Temporal discretization error: Proof of 

In order to show f J76|) . we first estimate 



E 



— X x 



B{X s )dW> 



A 



H 
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for all s,t G [0,T] with s < t and all K G N. More precisely, we have 



E 



* — X s 



B{X S ) dW t 



<5-E\\(e A ^ - I)X S \\ 2 H 



11 



+ 5-E 
+5-E 



A{t-u) 



F(X U ) du 



+ 5-E 



H 



l)B(X u )dWi 



K 



+5-E 



H 



e A ^B(X u )d(W u -W«) 
f{B{X u )-B{X.))dW* 

J s 



11 



H 



and using fP78|) shows 



E 



X t -X s - B(X s )dW< 



K 



H 



<5||MP (e A ^ - l)\\l (H) E\\(-Ay X s f H 



+ 5(t-s)Qf E\\e A{t ~ u) F{X u )\\ 2 H dt?J + 20i? 5 ^ sup 
+ 5^E||(e^-")-J)5(X u 



2 (l 



\HS(Uo,H) 



+ 5 (J\\\B(X U ) - B(X s )\\ 2 hs(u ,h) du 



for every s, t G [0, T] with s < t and every K G N. This shows 

2 



E 



B{X s )dW* 

<5R(t-sf J + 5{t-s) ( [ E\\F{X u )\\ 2 H du] + 20R 5 ( sup ^ 

\Js J \jeJ\J K J 



+ 5(£\\(-A)- 5 ( e A ^-l) 
+ 5R 2 E||X u -X,|&du) 



E 



HS(U ,H) 



du 
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and 



X t - X s - J B(X S ) dWi 

<5R{t-s) 21 + 5R{t-s) 2 + 20R 5 ( sup fiA 

\jej\j K J 



+ 5( / (t-u) 25 E (-A) d B(X 



HS(U ,H) 



du 



+ 5R 2 



-AY 



L(H) 



E\\X u -X a \\ 2 v du 



for every s,t E [0, T] with s < t and every K E N. Therefore, we obtain 



E 



Xt — X s 



B(X S ) dW r 



K 



<5R(t- s) 27 + 5R(t- sf 



H 



+ 20R 5 ( sup fjL^] + 5R ( I (t - u) 25 du] 
\jeJ\J K J \Js J 

+ 5i? 4 E\\X u -X s \\ 2 ^du^ 



and 
E 



ft 2 / \2o 

X t - X a - / 5(X S ) dW* < 10i? 3 (t - s) 27 + 20i? 5 ( sup fij ) 

# \jeJ\J K J 

+ 5R(t- s) {1+25) + 5i? 4 E ||X U - X S \\ 2 V ^ du^j 
for every s, t E [0, T] with s < t and every KeN. This shows 



E 



X t - X s - J B{X S ) dW« 



< 15R 3 (t-s 



1 2 7 



H 



+ 20R 5 I sup fiA 

\jeJ\J K J 



2a 




+ 5R 4 [ / [E\\X U -X. 



8\\vf,l dU 



and 
E 



Xt — X s 



B(X S ) dWi 



K 



( X 
< 15R 3 (t - s) 2j + 20R 5 sup ft 

H \jeJ\Jx J 

+ 5R* Qf (i2 (u - S )^^-^ 1 du ^ 
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and hence 



E 



Xt — X s 



B(X S ) dWi 



K 



< 15i? a (t - sf 1 + 20R b I sup ft 



H 



( SUp flj) 

\j&J\J K J 



2o 



+ 5R 5 



{u-s) 



min(2( 7 -/3),l) 



du 



for every s, t e [0, T] with s < t and every if £ N. Therefore, we obtain 



E 



X* — X x 



B(X S ) dW l u 



K 



H 



< 15R 3 (t - s) 27 + 20R 5 ( sup ft ) + 5R 5 (t 



( V a 

I sup Hj I 

\jeJ\J K J 



,(l+min(2( 7 -/3),l)) 



< 15i? 3 (t - s) 27 + 20i? 5 ( sup /i^ + 5i? 6 (t - s) 

\jeJ\J K J 



min(4(7-/3),2) 



and finally 



E 



Xt — X s 



B{X S ) dW t 



H 



<20i? 8 (t-s) min(4(7 - /3) ' 27) + 20i? 5 f sup fi,) (79) 

\jeJ\J K J 



for every s,t e [0, T] with s < t and every if e N. Now we prove (1761) . To 
this end note that 



E 



"i- 1 /.(M-l)h / ps \ 

V / e^^i?'^) [X s -X lh - B{X lh ) dW« dWf 



(=0 

< V / E £?'(X ih ) X s - X i?i - / S(Xjh) dW* 

i=0 V .//ft / 



HS(U ,H) 



, m—\ 



< 



1=0 



(l+l)h 



E 



K 



Ih 



ds 



H 
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holds for every m G {0, 1, . . . , M} and every M, if el Hence, (lT§j) yields 

m -i f(l+l)h f rs \ 

/ e^-^i^X,,) [X s -X lh - B(X lh ) dW* dW t 

l=Q Jlh V </«/i / 

< 20/? 10 V / ( S - ihf^-^) + ( sup ^ ) ds 



(m— 1 „ 
y / 
z=o 



m -! Al+l)h \ / \2a 

" +20iT ao T( sup ^ 

/ \jeJ\J K 



and therefore 

m-l 



m-i n(l+l)h / rs \ 

V / e^—^S'^) [X s -X lh - B{X lh ) dW« dW ( 

l=Q Jlh V Jlh J 

<20R 10 Mh( 1+rain ^-^> 2 ^ + 20R U ( sup ^Y" 

V jeJ\J K J 

< 20 J R 11 /i min(4(7 - /3) ' 27) + 20 J R 11 f sup n}\ 

\jej\j K ) 



ii 



and finally 

m -! »(J+l)fc 



E 



V / e^^B'^) [X s -X lh - B{X lh ) dW« dW ( K 

l=Q Jlh V ./Z/i / 



20i? 13 n / 

^ Mmi n(4( 7 - /3 ), 27 ) + 20i? ( SU P 



j&J\J K 



for every m G {0, 1, . . . , M} and every M, K G N. 
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6.6 Lipschitz estimates: Proof of ( 1771 ) 

Before we estimate E \\Z% M ' K - Y^' M ' K f R for m G {0,1, ... ,M} and for 
N, M, K e N, we need some preparations. More precisely, we have 



E 



B'(X lh ) f B(X lh ) dW* - B'(Y^ ui< 

Jlh 

( 



Ih 



B[Y l w )dW K 



HS(U ,H) 



E 



B'(X, 



ih) 



v r B(x lh ) gj d( gj ,w u ) 

tTZ Jlh 



u 



/ 



^/ lyN,M,K 



N.M.K 



9jd(gj,W u ) 



u 



HS(U ,H) 



and hence 
E 



S'(Jf Jfc ) I* B{X lh ) dW« - B' (y^*) f b(y; 
Jih v ; Jih K 

( 



HS(U ,H) 



E 



//J 



\ 



— b'(y n,m,k 



\ 



Y,B{Y l N ^ K ) 9] {^W s -Wih) l 



HS(U ,H) 



and 
E 



Jih v 



Ih 



B[Y^ M ' K )dW K 



HS(U ,H) 



E 



jeJ K 



u 



HS(U ,H) 
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for every s E [lh,(l + l)h], I E {0, 1, . . . , M - 1} and every M,K E N. 
Moreover, we have 



E 



B\X lh ) £ B(X lh ) dW u K - B' (Y/w) £ B (Y^) dW u K 



HS(U ,H) 



B'(X lh ) (B(X lh ) 9j ) 



b'(y^ k ) (b(y^ k ) 9j ) } ( 9j , w s - w lh)u 



HS(U ,H) 



and 



E 



\X lh ) f B{X lh )dW^-B'(Y l N ^ K ) f B{Y t N ^ K )dW^ 

J lh J lh 

= y £E\\B'(X lh )(B(X lh )g s )-B , (Y l 



HS(Uo,H) 



rN,M,K 



B[Y, | gj 



HS(Uo,H) 



E\(g j: W s -W lh ) l 



for every s E [lh, (I + l)h], I E {0, 1, . . . , M - 1} and every M, K E N. This 
shows 



B'{X lh ) £ B(X lh ) dW« - B' (¥»>">*) £ B (¥»>">*) dW« 
= £ ^ r E\\B\X lh mX lh ) 9j )-B'(Y^ M ' K )(B(Y^ K ) 9] 



and therefore 



HS(U ,H) 
2 

HS{U ,H) 

s-lh) 
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< ^^\B\X Ul ){B{X lh )g 3 )g k -B'(Y^ K )^ 



j,keJ 



HS(U ,H) 
2 



H 



(s - lh) 
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for every s E [Ih, (I + i)h], I E {0, 1, . . . , M- 1} and every M, K E N. Hence, 
we obtain 



E 



'N.M.K 



B'{X lh ) [ S B(X lh ) dW« - B'(Y^ 

Jlh V 

= (s -lh)- E \\B'(X lh ) B(X lh ) - B'(Y^ M ' K ) B(Y^ M ' K ) 
and finally 
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< R -E 



Xm-Yi 



N,M.K 



HS(U ,H) 

\ < 80 > 



for every s E [lh,(l + l)h], I E {0,1,..., M — 1} and every M,K E N. 
Additionally, ( |69|) shows 
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and hence 

£ \\^N,M,K _ yN,M,K\ 
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for every m e {0, 1, . . . , M} and every N, M, K EN. Moreover, we have 
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HS(U ,H) 



and due to (1HU1) we obtain 
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for every m G {0, 1, . . . , M} and every N, M, if el. Finally, we obtain 



'm—l 



j£ \ \r?N,M,K _ y n M,K\\' 2 



H <9R 3 h E E \\X lh -Y l 



N,M,K 



9R 



1=0 

'm—l 



N,M,K 



, 1=0 



for every m G {0, 1, . . . , M} and every N, M, K G N. 



6.7 Iterated integral identity: Proof of 

First of all, we have 



(m + l)T 
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iiiT 
A/ 
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P-a.s. and hence 
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mT I mT 

M M 



d(g k ,W u ) u d(g j ,W s ) l 



P-a.s. for all m G {0, 1, . . . , M — 1} and all N,M,K G N. Moreover, since 
the bilinear operator B' (Y£ > m ' k ) B (Y*' m > k ) G H (U q , H) is symmetric 
(see Assumption [3]) and since 



(m + l)T 
Af 



mT / mT 

A/ M 



.2 TUj 



P-a.s. holds for all j G J K , m G {0, 1, . . . , M - 1} and all N, M, K G N (see 
(3.6) in Section 10.3 in [36J), we obtain 
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P-a.s. for all m G {0, 1, . . . , M - 1} and all N,M,K G N. The fact 
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M M M M 
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for all j G J K , m G {0, 1, . . . , M - 1} and all N,M,K G N (see (3.15) in 
Section 10.3 in [36]) then yields 



r - B\Y»> M > K )( [ B(Y»^ K )dwAdW? 

mT \ I mT / 

1 B^^^ 



2 



T 

~ 2M 



for all m G {0, 1, . . . , M - 1} and all N, M, K G N which shows ([65]). 
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